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ABSTRACT 

The  application  of  error  detection  and  correction  codes 
has  advanced  markedly  with  the  advent  of  digital  technology. 
However,  the  strides  made  towards  employing  the  techniques, 
encoding  and  decoding,  properly  have  been  rather  limited 
by  the  Naval  forces  in  the  United  States.   This  paper  deve- 
lops a  computer  program,  which   if  utilized  properly,  would 
aid  in  deciding  what  error  correcting  scheme  is  best  suited 
for  a  specific  channel. 

The  results  obtained  from  testing  a  rate  1/n  convolu- 
tional  code,  over  a  simulated  channel,  using  a  Viterbi 
decoder  shows  that  this  is  an  effective  analysis  procedure. 
Though  the  test  runs  were  lengthy,  much  of  the  time  required 
was  for  noise  simulation.   This  would  not  be  a  factor  if  actual 
channel  noise  recordings  had  been  available. 
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I.   INTRODUCTION 

The  distortion  of  a  digital  transmission,  caused  by 
noise,  may  require  that  techniques  be  used  allowing  for  the 
detection  and  eventual  correction  of  errors.   In  digital 
communications,  the  most  significant  performance  parameter 
from  either  a  bit  or  message  standpoint  is  that  of  probabi- 
lity of  error.   A  means  of  minimizing  the  probability  of 
error  is  therefore  essential  to  effective  communication. 

Error  detection  and  correction  was  given  great  impetus 
by  Shannon's  paper  of  1948,  [Ref . 7  ],  which  extended  the 
promise  of  reliable  recovery  of  digital  data  perturbed  by 
noise  (Shannon's  Second  Theorem).   The  noisy  coding  theorem 
provides  that  messages  can  be  transmitted  with  arbitrarily 
small  error  if  the  source  rate  is  no  more  than  the  channel 
capacity.1   While  this  sets  a  goal  and  the  conditions  neces- 
sary to  attain  it,  the  precise  method  for  obtaining  a  speci- 
fic scheme  was  not  set  forth.   However,  Shannon  did  use  a 
random  coding  scheme  in  proving  this  theorem.   As  was  shown 
by  Rice  in  1950,  [Ref.  6  ],  following  Shannon's  example, 
choosing  codewords  randomly  leads  to  the  result  that  as  the 
code  becomes  very  long,  channel  capacity  is  approached  and 


1  Channel  capacity  is  the  maximum  average  information 
input  to  a  channel,  which  means  it  is  properly  matched  to 
the  channel  less  the  average  uncertainty  at  the  receiver 
due  to  channel  noise. 
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the  probability  of  error  can  be  made  to  decrease  exponen- 
tially with  code  length.  But  choosing  codewords  randomly 
is  not  a  practical  scheme. 

Until  the  1960's  the  application  of  channel  coding  theory 
was  slow  in  its  development.   The  establishment  of  digital 
circuit  technology  and  the  realization  of  a  theoretical 
channel  (satellite  communication;  AWGN)  provided  a  powerful 
stimulus  for  the  utilization  of  practical  error  detection  and 
correction  techniques.   The  technological  improvements  have 
continued  until  today  fairly  large  memories  are  held  by 
small  chips,  and  the  components  needed  to  implement  powerful 
coding  schemes  are  available  in  reduced  sizes  and  at  reason- 
able prices. 

The  objective  of  this  paper  is  to  demonstrate  the  bene- 
fits derived  from  a  minicomputer  (DEC  PDP-11/40)  simulation 
of  channel  noise  applied  to  a  specific  error  detection  and 
correction  scheme.   This  simulation  allows  for  the  determina- 
tion of  a  probability  of  bit  error  for  a  specifiable  noise 
density.   Results  from  the  simulation  are  then  used  to  deter- 
mine statistical  trends  of  a  particular  error  correction 
technique  for  various  internal  coding  parameters. 

The  error  detection  and  correction  scheme  chosen  by  the 
author  is  the  Viterbi  decoding,  [Ref.3  ],  algorithm  on  a  rate 
1/n  convolutional  code.   Both  the  encoding  and  decoding 
technique  are  of  current  interest  in  extraterrestrial  com- 
munications, with  the  Viterbi  algorithm  being  principally 
considered  for  satellite  communications.   Therefore,  in 
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order  to  provide  a  complete  presentation,  the  scope  of  this 
thesis  will  encompass  the  basic   principles  of  convolutional 
encoding  and  Viterbi  decoding,  as  well  as  the  structure  of 
the  simulation  program,  results  and  conclusions. 

All  programming  was  accomplished  on  the  DEC  PDP-11/40 
in  machine  language,  thus  presenting  the  opportunity  to 
gain  valuable  insight  into  the  computer  structure  while 
completing  this  thesis. 
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II.   BACKGROUND 

Digital  communication  systems  (figure  1)  are  usually 
designed  to  minimize  the  probability  of  error  of  the  re- 
ceived data  (bits) ,  introduced  by  noise  (gaussian,  burst, 
fading,  multipath,  etc.).   Thus,  error  detection  and  correc- 
tion may  be  varied  according  to  the  degree  and  type  of  noise. 
There  are  two  common  forms  of  error  protection:  (1)  retrans- 
mission and  (2)  controlled  redundancy.   The  first  may  be 
applied  to  systems  which  are  not  critically  affected  by  the 
ensuing  delay,  whereas  the  second  does  not  require  signifi- 
cant delays.   This  thesis  deals  with  the  latter  in  an 
attempt  to  gain  statistical  knowledge  of  the  particular 
coding  scheme  (controlled  redundancy)  under  study. 

Controlled  redundancy  techniques  are  commonly  divided 
into  two  groups:  (1)  block  codes  and  (2)  convolutional 
(tree)  codes.   Again,  owing  to  the  characteristics  of  the 
channel,  one  of  these  two  codes  may  be  chosen,  along  with 
a  suitable  decoding  scheme,  to  achieve  the  desired  probabi- 
lity of  error.   A  decoding  scheme  for  block  codes  is  batch 
oriented,  as  is  its  encoding,  using  definite  algebraic 
operations  related  to  the  segmented  structure  of  the  code. 
On  the  other  hand,  convolutional  codes  are  decoded  by  a 
statistical  procedure  due  to  its  continuous  (bit  by  bit) 
nature. 
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A.   PRINCIPLES  OF  CONVOLUTIONAL  ENCODING 

In  1955,  P.  Elias  first  proposed  the  use  of  convolutional 
(tree)  codes  for  the  discrete  time  memoryless  channel, 
[Ref  2]  .   This  technique  extended  the  promise  of  providing  a 
class  of  codes  (linear)  whose  performance  would  prove  supe- 
rior to  that  of  block  codes  of  the  same  length,  [Ref. 9  ] . 
The  development  of  these  codes  also  gave  promise  of  provid- 
ing for  a  decoder  complexity  increasing  no  more  than  linearly 
with  block  length  and/or  encoder  memory.   These  conjectures 
have  been  verified,  for  the  most  part,  by  such  contributors 
as  Fano,  Reiffen,  Forney  and  Berlekamp.  [Ref.  4] 

Since  convolutional  codes  form  a  definite  discipline 
in  coding  theory,  there  are  some  aspects  that  will  require 
clarification  in  conjunction  with  the  implementation  of  the 
code.   The  next  few  subsections  will  accomplish  this  task. 

1.   Description 

The  process  of  block  encoding  segments  (blocks)  of 
an  indefinite  length  input  sequence  into  code  blocks  pro- 
ceeds in  a  discrete  manner,  i.e.,  the  contents  of  one  block 
has  no  bearing  on  the  encoding  of  another  block.   The  code 
blocks  may  be  completely  changed  from  the  input  block  or 
they  may  include  the  input  block  with  check  bits  determined 
by  a  definite  algorithm.   However,  when  an  input  sequence 
is  processed   through  a  convolutional  encoder  the  input/ 
output  is  continuous  (sequential  bit  to  bit  dependence  in 
the  encoder)  and  the  output  is  generated  at  a  specifiable 
number  of  code  bits  for  each  encoder  input  bit. 
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The  term  "convolutional"  originates  from  the  observation 
that  the  encoded  sequence  can  be  regarded  as  a  convolution2 
of  the  input  sequence  with  related  generator  sequences.   An 
example  of  such  a  procedure  is  the  binary  convolution 
(denoted  by  *)  of  two  sequences,  x  and  g_,  where  (+)  and  •  are 
binary  addition  and  multiplication,  respectively. 
Input  sequence: 

x(tQ)  =  (xQx1x2x3x4  . .  .)  =  (1  0  1  1  0  .  .  .)  ; 

Generator  sequence: 

£(tQ)  =  (gQg1g2g3g4  ...)  =(10011...), 

where  all  bits  in  x   and  g_  are  zero  for  x  .   and  g  .  . 
Then  the  output  (code)  bit  formed  by  the  convolution  of 
x(tQ)   with  g(tQ)   is 

Z(tQ)  =  x(t0)*  2(t0)  =  0  . 
For  time   t„  +  At,   the  sequences  and  output  are 

x(tQ+At)  =  (xQx;Lx2x3x4  ...)  =  (1  0  1  1  0  ...), 

g_(tQ+At)  =  (gog;Lg2g3g4  ...)  =  (1  0  0  1  1  ...), 

y(tQ+At)  =  (xQ  •  gQ)  =  (1  •  1)  =  1  . 
For  time   t_  +  2At,   the  sequences  and  output  are 


2The  convolution  of  two  functions,  f,  and  f„,  is  equal 
to  /■» 

f1(t-T)f2(T)dT  , 
/  —oo 

which  means  one  function  is  folded  in  time,  then  their 
product  is  integrated  as  one  is  shifted  pass  the  other. 
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x(tQ+2At)    =    (XgX^   x   x      . . .)    =    (1    0    1    1    0    . . .) , 

2(tQ+2At)    =    (gog;Lg2g3Q4   .  ... }    =    (1    0    0    1    1    ...), 

y(t0+2At)    =    (xQ    ■    gx)©  (x1    •    gQ)    =0. 

And  the  same  procedure  is  applied  to  subsequent  increments 
of   At  (one  bit  shift) . 

This  technique  extends  the  concept  of  block  encoding 
to  permit  memory  from  bit  to  bit  continuously  as  against 
memory  within  the  block.   Using  the  operation  shown  above, 
a  conventional  convolutional  encoder  (figure  2)  may  be 
defined  as  a  linear  sequential  machine  with  k-inputs  and 
n-outputs,  where   n  >  k,  usually.   This  machine  is  constant, 
linear,  causal,  and  finite  state  with  operations  over  a 
finite  field   F,   commonly  binary. 

The  parameter  k   is  the  number  of  input  bits  which 
enter  the  encoder  in  a  time  increment   (At) .   These  bits 
along  with  others  retained  in  the  encoder  memory  (finite 
length) ,  form  n-output  bits  (code  bits) .   For  the  remainder 
of  this  paper   k   will  have  a  value  of   1,   but  it  should 
be  remembered  that  this  parameter  can  take  on  other  values, 
as  can  be  seen  in  figure  3. 

The  definition  specified  a  finite  state  machine,  where- 
as the  example  dealt  with  indefinite  length  sequences. 
Since  an  infinite  length  memory  is  not  practical  it  becomes 
necessary  to  limit  the  basic  concept  (generator  sequence) 
and  define  some  new  expressions  employed  in  the  actual 
encoder  configuration. 
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Figure  2.   Conventional  (n,k)  Convolutional  Encoder 
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a.  Code  Rate  (k/n) 

The  code  rate  is  an  expression  of  the  number 
of  k-input  (1)  bits  per  n-output  bits,  that  is,  k/n  (1/n) . 
This  value  may  be  considered  a  very  close  approximation 
to  the  actual  value.   The  actual  value  would  also  include 
a  very  small  number  of  zero  code  bits  used  to  terminate  the 
code  sequence. 

b.  Memory  Length  (M) 

For  a  practical  encoder  £  must  be  limited  in 
length.   The  memory  length  is  a  measure  of  this  length  and 
is  the  minimum  number  of  memory  cells  required  to  generate 
a  code.   A  representation  of  an  encoder,  with  a  memory 
length  of  2,  is  shown  in  figure  4.   This  configuration 
corresponds  to  the  Mealy  machine  in  automata  theory. 

c.  Constraint  Length  (K) 

A  more  conventional  expression  of  the  length 
of  time  which  an  input  bit  affects  the  output  sequence  is 
that  of  constraint  length.   This  term  is  simply  the  memory 
length  plus  one,  and  is  depicted  as  the  number  of  memory 
elements  in  a  Moore  machine  configuration  of  an  encoder 
(figure  5) .   As  may  be  expected  this  value  is  one  factor 
in  determining  the  complexity  of  the  error  detection  and 
correction  scheme. 

d.  Free  Distance  (df) 

The  term  distance,  as  applied  to  coding,  refers 
to  the  number  of  differing  bits  between  two  code  sequences 
of  the  same  length,  as  seen  below: 
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code  sequence  1:  (1001011001) 
code  sequence  2:  (1011001001) 
Differences  between  the  two  sequences  exist  in  the  third 
and  sixth  bits  only,  therefore  the  distance  between  the  two 
sequences  is  2.   When  using  block  codes  to  construct  a  code 
sequence  the  minimum  distance  between  all  pairs  of  codewords 
(code  blocks)  is  a  definite  indication  of  the  error  detect- 
ing and  correcting  capability  of  the  code.   However,  if  the 
code  does  not  employ  a  block  configuration,  the  use  of  this 
value  may  not  be  justified  to  indicate  the  codes  error 
detecting  and  correcting  capabilities. 

This  is  the  case  when  discussing  convolutional  codes. 
Free  distance  is  the  term  for  the  value  expressing  the 
theoretical  error  correcting  capability  of  a  convolutional 
code.   This  value  is  defined  as  that  minimum  number  of  set 
bits  occurring  in  a  code  sequence  which  resulted  from  the 
input  of  a  nonzero  sequence.   Table  I  shows  how  the  free 
distance   (df)  can  be  determined  for  a  rate  1/2  convolu- 
tional code  with   K=3.   As  can  be  seen,  the  larger  the 
value  of   K   then  the  number  of  possible  input  sequence  to 
be  considered  also  increases  at  a  great  rate.   But  with  the 
aid  of  computers  the  tedious  grinding  process  is  accomplished 
rather  quickly. 

The  convolutional  encoder  is  constructed  to  take  full 
advantage  of  the  free  distance  applied  to  a  specific  code 
rate  and  constraint  length.  This  is  accomplished  through 
the  proper  selection  of  generator  sequences.   However,  it 
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should  be  noted  that  with  the  computers  currently  available 
the  drudgery  can  be  passed  on  to  the  machine. 
2.   Properties 

The  properties  of  the  encoder,  as  related  to  the 
generator  sequence,  are  discussed  in  this  section.   For 
convenience,  the  set  of  generator  sequences   (g_,  ,  .  .  .  ,g_  ) 
used  to  generate  the  convolutional  code  is  denoted  by 
[G] ,  a  matrix  representation  of  the  encoder  discussed  in 
the  next  section.   Figure  6  is  useful  in  relating  the 
following  discussion  to  a  physical  communication  system. 

a.  Property  1 

Foremost,  the  useful  encoder  should  generate 
a  code  which  will  yield  the  fewest  number  of  errors  in  the 
codeword  estimator.   The  estimator  is  a  demodulation  scheme 
at  the  receiver  determining  as  accurately  as  possible  the 
received  sequence  before  actual  decoding  takes  place.   Along 
with  this  error  minimization,  the  complexity  of  the  estima- 
tor requires  minimization  to  develop  a  useful  error  correc- 
tion system.   Error  minimization  and  complexity  minimization 
together  may  require  some  compromise  in  the  practical 
system  depending  on  the  desired  reception  quality  and  sys- 
tem environment. 

b.  Property  2 

The  encoder  [G]  must  have  an  inverse  relation 
for  decoding  purposes.   This  is  realized  by  the  relation 

yjG]"1  =  x[G]  [G]"1  -  DPx 
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for  all   x.   The  matrix  (decoder),  [G]~    is  physically 
realizable  and  is  a  pseudo  inverse,  in  that,  the  expression 
[G] [G]   ,   yields  a  delay   Dp,   where   p   is  the  number 
of  time  intervals  before  decoding  occurs  in  the  expression. 
If  p   is  zero  the  decoding  occurs  on  receipt  of  the 
sequence  and  the  above  relation  becomes 

yjGj"1  =  x[G]  [G]"1  =  x, 

where   [G]    =  [G]     when   p  =  0. 

A  mention  of  the  possibility  of  catastrophic  error 
propagation,  [Ref   ] ,  is  required  at  this  time  to  stress 
the  point  that   [G]    must  be  feedback-free.   This  may  be 
interpreted  as  meaning  that  the  n   generator  sequences, 
in  polynomial  form  should  not  have  a  common  factor.  If  a 
common   factor  exists  then  the  input  of  a  sequence  with  a 
finite  number  of  set  bits  to  the  encoder,  may  be  decoded, 
after  noise  is  added,  as  a  sequence  with  an.  infinite 
number  of  set  bits.   This  is  catastrophic  error  propagation 
and  is  avoided  by  making  [G] ,  thus  [G]    feedback-free. 
c .   Property  3 

The  matrix   [G]   must  also  meet  the  requirements 
mentioned  in  the  definition  of  a  convolutional  encoder. 

(1)  Constant  (time  invariant) .   The  time 
invariance  of  the  encoder  is  represented  by 

[G]  (DPx)  -  DPy_  , 
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which  means  if  all  the  inputs  are  shifted  in  time,  then 
all  the  outputs  are  shifted  accordingly. 

(2)  Linear.   The  output  sequence  resulting 
from  the  superposition  of  two  input  sequences,  x,      must 
be  equal  to  the  superposition  of  the  two  output  (code) 
sequences,  y_,  that  would  result  from  the  inputs  entered 
separately.   This  is  also  necessary  for  the  multiplication 
of  an  input  by  a  scalar  as  shown  below.   If   G:x-»-y_,   then 
G(x1+x2)  =  GCx^)  +  G(x2)  =   Y-i   +   Y.2'       and   G(ax-,)  =  aGfx,)  = 
ay-,   where   a   is  an  element  in  the  field   F  (usually  GF(2)) 

(3)  Causal.   The  existence  of  a  nonzero  output, 
y,   prior  to  the  input  of  a  nonzero  sequence,   x,   is  for- 
bidden.  This  is  accomplished  by  the  encoder  output  being 
zero  when  its  memory  elements  are  all  zero,  therefore  zero 
inputs. 

(4)  Finite  State.   The  states  of  the  encoder 

are  finite  in  number  due  to  the  value  of  the  constraint 

K— 1 
length.   Thus,  a  binary  encoder  has   2      possible  states. 

Each  state  being  those  input  elements  involved  in  the 

generation  of  the  next  output  (code)  bit,  when  combined 

with  the  next  input  bit  entering  the  encoder. 

d.  Property  4 

The  chosen  encoder  [G]  should  have  the  minimum 
number  of  memory  elements  to  generate  the  code. 

e.  Property  5 

The  chosen  encoder   [G]   should  generate  a  code 
yielding  the  fewest  number  of  decoding  errors  per  error 
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event  (channel  noise) .   Thereby,  the  fastest  correct 
decoding  decision  is  made  at  the  receiver. 

The  preceding  paragraphs  provide  a  basis  by  which  the 
acceptable  class  of  encoders  is  greatly  reduced  in  number. 
However,  the  selection  of   [G]   in  its  optimum  form  for  a 
given  channel  is  a  long  drawn  out  process,  possibly  best 
suited  for  computer  analysis.   Theoretical  selections  have 
been  made  by  many  people,  [Ref.   ],  who  based  their  deci- 
sions on  the  free  distance  alone. 

3 .   Representation 

The  encoder  configurations  used  to  this  point  are 
accurate  physical  descriptions,  but  they  lack  a  convenient 
form  needed  for  analysis.   There  are  five  encoder  repre- 
sentations which  show  varying  degrees  of  the  code's  struc- 
ture.  Each  is  discussed  in  the  following  sections, 
a.   Polynomial 

Using  the  D-transform,  the  input  and  output 
sequences  appear  as  polynomials  of  degree  much  greater  then 
the  constraint  length   (K) . 

-1  2 

Input,   X(D)  =  ...  +  x_1D    +  xQ  +  x^p    +    x2D   +  ... 

output,  Yi(D)  =  ...  +  y±  jp'1  +  yi/0  +  YitiD 


.2 
1,2 


+  y .   D   +  .  .  .  , 


where   1  £  i  &   n.   The  generator  sequence,  for  the  Mealy 
machine  configuration  in  figure  7,  is  a  polynomial  of 
degree  equal  to   K-l.   The  mapping  of   X  ■*  Y    is  accom- 
plished by  polynomial  multiplication. 
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Y.  (D)  =  I     (x.  •  g.    ,)D! 


where   L   is  the  bit  length  of  the  input  sequence.   This 
representation  has  little  advantage  in  a  code  structure 
analysis,  but  it  is  very  definitive  as  to  the  procedure  for 
code  generation. 

b.   Matrix 

A  matrix,  [G]   representation  of  a  rate  1/n 
encoder  starting  at  state  zero  and  time  zero  is  achieved 
by  the  manipulation  of  each  generator  sequence.   When 
y  =  [G] x  , [G]   may  be  shown  as  a  matrix  whose  row  vectors 
are  the  generator  sequences  (figure  8a) .   Using  this  nota- 
tion restricts   x   to  be  a  length   K.   However,  the  require- 
ment to  change   x   for  each  code  computation  is  eliminated 
by  the  matrix  in  figure  8b.   This  notation  incorporates 
the  shifting  of   x  past  each  2-  d  *•  i  £   n)   into  the 
matrix   [G] .   For  example,   (figure  8  follows  the  example), 


31 


X  = 


1 

0 
0 

1 
0 

1 
1 
1 

0 

1 
1 

1 

0 
0 

1 

0 

0 


n  =  2 
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2j=  (l  0  l) 

22=  (1  1  1) 
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This  form  allows  for  a  fast  generation  of  the  code  and 
is  definitely  a  model  which  could  be  implemented  on  a  com- 
puter with  little  difficulty.   Again,  there  is  a  lack  of 
insight  into  the  code  structure  except  for  the  possible 
application  of  matrix  algebra  theory. 
c.   Tree 

The  most  common  representation  is  a  tree  dia- 
gram (figure  9)  which  incorporates  the  branch  and  nodal 
properties  of  the  code  generation.   The  base  branch  corres- 
ponds to  the  initial  state  of  the  encoder  prior  to  a  non- 
zero entry.   The  first  node  (a)  is  state  zero  for  the  en- 
coder.  When  a   1   begins  the  sequence  the  first  lower 
branch  is  chosen  and  the  code  is  found  on  this  branch 
leading  to  node  (b)  or  state   1.   The  subsequent  input 
bits  dictate  whether  the  tree  is  followed  up  or  down, 
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Figure  8.   Generalized  Matrix  Representation  of  a  Rate 
1/2  Convolutional  Encoder   (K=3) . 
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corresponding  to  inputs  of   0   and   1   respectively.  If 
this  is  continued  to  the   K    branch,  all  of  the  possible 
input  sequences  of   K  bits  have  been  represented,  thus 
each  possible  encoder  state.   Beyond  the   K    branch  it  is 
readily  recognized  that  the  code  symbols  on  the  branches 
leaving  from  the  two  nodes  labeled  ( a)  are  identical  iden- 
tical to  those  at  the  base  of  the  tree.   Even  though  the 
structure  has  become  recurrent,  for  long  sequences  this 
diagram  may  become  too  large  to  analyze  easily. 
d.   Trellis 

Since  the  input  sequence   (1  0  0  ...)  and 
(0  0  0  ...)   generate  the  same  n-bit  codewords  after  the 
third  branch  (figure  9) ,  then  both  nodes  labeled   a   can  be 
joined  together.   The  recursive  character  of  the  tree 
diagram  lends  itself  to  be  redrawn  with  remerging  branches, 
thus  forming  a  trellis  (criss-cross)  diagram  (figure  10) . 
The  encoder  state  is  the  basis  for  the  trellis  diagram. 
These  states  represent  those  input  bits  in  the  memory  cells 

which  will  generate  the  code  bits  when  the  next  input  bit 

K— 1 
enters  the  encoder.   Therefore,  there  are   2      states  on 

either  side  of  the  diagram  representing  the  present  and 

next  states.   The  branches  connecting  the  states  denote  the 

code  bits  being  generated  when   a   1   (dashed  line)  or   a 

0   (solid  line)  enters  the  encoder.   Again,  after  the  third 

input  bit  the  complete  trellis  diagram  is  specified  and  is 

repeated  for  all  succeeding  branches. 
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Figure  9.   Tree  Representation  of  Convolutional 
Encoder  in  Figure  5. 
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For  short  constraint  lengths   (K  <  10) ,   this  repre- 
sentation is  highly  desirable  due  to  its  compactness.   The 
trellis  diagram  is  used  in  the  Viterbi  decoding  algorithm 
and  easily  lends  itself  to  be  listed  in  a  computer  program. 
e.   State  Diagram 

The  last  encoder  representation  is  similar  to 
the  trellis,  but  is  spread  out  more  to  allow  more  parameters 
to  be  placed  on  the  branches  between  states  (nodes) .   The 
state  diagram  (figure  11)  has  the  input   0   and   1   denoted 
in  the  same  manner  as  for  the  trellis,  but  some  additional 
notation  is  placed  on' the  branches. 

(1)  D^.  The  value  of  q  denotes  the  weight 
of  the  code  bits  for  that  branch. 

(2)  L^.  The  value  of  r  denotes  the  branch 
length  of  a  path  from  state   x   to  state   y. 

(3)  N  .  The  value  of  s  denotes  the  number 
of  input  one  branches  encountered  in  a  path  from  state  x 
to  state   y. 

All  of  the  above  notation  is  used  in  figure  11  and  the 
results  of  a  state  path  for  the  input  sequence   (10  10  0) 
are  shown. 

Since  the  state  diagram  is  a  directed  graph,  a  transfer 
function  can  be  determined  using  the  theory  of  directed 
graphs,  [Ref.  4,  pp  2A-1-2A-11] .   The  transfer  function  con- 
sists of  various  powers  of  the  three  measurements  listed 
above.   These  values  are  used  to  determine  the  properties 
(D,  L,  N)  for  all  paths  in  the  convolutional  code,  when  the 
transfer  function  is  in  expanded  form. 
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Solid  Line  -  Input  0 
Dashed  Line  -  Input  1 


DLN 
X  =  (10100);  Path  Parameters  =  DLN 


Figure  11.   State  Diagram  Representation  of 

Convolutional  Encoder  in  Figure  5. 
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The  trellis  and  state  diagrams  are  the  most  complete 
representations  and  offer  the  analyst  a  chance  to  study 
the  path  structure  of  the  convolutional  code.   Each  diagram 
yield  a  compact  representation,  easily  arranged  for  small 
values  of   K.   However,  the  looping  branches  necessary  for 
the  construction  of  a  state  diagram  become  very  confusing 
and  difficult  to  arrange  for  larger  constraint  lengths. 
Therefore,  the  trellis  diagram  is  not  only  easy  to  construct 
and  understand,  but  if  the  values  of   D,  L,  N,   are  needed 
for  analysis  they  can  incorporated  into  the  diagram  in  an 
orderly  fashion. 

The  discussion  of  convolutional  encoding  principles 
was  basic  and  is  by  no  means  a  complete  detailed  course  in 
convolutional  codes.   However,  the  purpose  of  this  section 
is  to  introduce  the  reader  to  concepts  which  will  aid  him 
in  his  understanding  of  the  thesis  computer  program. 

B.   PRINCIPLES  OF  VITERBI  DECODING 

Decoding  is  the  inverse  operation  to  encoding  and  is 
intended  to  recover  the  source  bits  with  all,  or  almost 
all,  of  the  channel  errors  removed.   The  decoder,  which 
may  be  implemented  as  hardware  or  software,  utilizes  the 
encoded  bits  to  detect  and/or  correct  errors.   Error  detec- 
tion is  similar  in  complexity  to  the  encoding  operation. 
Error  correction,  however,  is  inevitably  a  more  complicated 
process  than  encoding,  since  the  goal  is  to  reduce  the 
probability  of  a  decoded  bit  error. 
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A  scheme  for  decoding  convolutional  codes  was  proposed 
by  Viterbi  in  1967,  [Ref.  8].   Trt  was  shown  by  Viterbi  that 
an  optimum  decoding  procedure  existed  for  a  statistically 
independent  (bit  to  bit)  input  sequence  transmitted  over  a 
channel  whose  errors  occur  independently  from  channel  bit 
to  channel  bit.   Subsequently,  Forney  found  that  the  Viterbi 
decoding  algorithm  is  synonomous  with  maximum  likelihood 
sequence  decoding,  [Ref.  3]. 

1.   Description 

Simply  stated,  the  Viterbi  algorithm  is  a  solution 
to  the  problem  of  finding  the  most  likely  encoded  sequence 
through  a  state  (finite)  diagram  representing  the  encoder. 
For  a  decoder  to  minimize  the  overall  error  probability 
of  a  decoded  bit  by  brute  force,  maximum  likelihood  decoding 
would  mean  calculating  the  likelihood  of  the  received 
sequence  on  all  paths  of  the  encoder  state  diagram.   How- 
ever, there  are  two  factors  which  reduce  the  complexity  of 
this  problem.   The  first  is  the  fixed  periodic  structure 
of  the  encoder  trellis  diagram,  and  the  second  is  the  code 
characteristic  of  remerging  paths  after  the  same   K   input 
bits  are  applied  to  two  different  paths  (figure  12) .   For 
these  reasons,  the  trellis  diagram  is  an  ideal  tabular 
representation  of  the  flow  of  the  code  at  any  instant. 

The  decoding  process,  as  mentioned  above,  was  found 
optimal  for  a  statistically  independent  input  sequence  in 
discrete  time.   Along  with  this  stipulation,  the  channel 
noise  is  also  memoryless,  as  is  the  case  for  a  binary 
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symmetric  channel  (BSC)  with  hard  decision  demodulation 
(2   levels)  or  an  additive  white  aaussian  noise  (AWGN) 
channel  with  soft  decision  demodulation  (2   levels,  Q  >  1) . 
The  BSC  errors  transform  a   0   to  a   1   and  a   1   to  a   0 
and  occur  independently  from  bit  to  bit  with  probability 
p.   In  the  AWGN  channel,  the  probability  of  a  given  quan- 
tized value  (0  to  2-1)  of  a  received  bit  is  determined 
from  the  gaussian  probability  density  function.   These 
values  of  the  received  bits  are  used  to  determine  the  most 
likely  received  sequence. 

A  scoring  procedure  is  employed  indicating  the  trellis 
path  which  shows  the  least  difference  from  the  received 
sequence.   This  is  accomplished  by  using  the  concept  of 
distance  for  the  BSC  and  a  numerical  difference  for  the 
AWGN  channel.   The  latter  may  be  called  the  innerproduct 
of  the  received  bit  and  the  calculated  (trellis)  trans- 
mitted bit.   This  is  implemented  digitally  by  taking  the 
exclusive  OR  of  the  two  Q  bit  representation  of  the  bits 
above.   Both  (BSC  and  AWGN)  techniques  yield  a  value  that 
can  be  used  to  determine  the  most  likely  (lowest  score) 
path. 

2 .   Implementation 

In  order  to  describe  effectively  the  implementation 
of  the  Viterbi  decoding  algorithm  an  example  follows  with 
step  by  step  explanations  supplemented  by  appropriate  dia- 
grams.  This  example  will  be  for  a  rate  1/2  convolutional 
code  with   K  =  3.   The  generator  sequences  are  5g (1  0  1) 
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and   7R  (1  1  1) .   The  following  sequence  is  the  encoder  in- 
put:  (101100011Q100010  ...}. 

The  corresponding  code  sequence  is 
(11  01  00  10  10  11  00  11  10  10  00  01  11  00  11  01  11  00  ...). 
An  error  sequence  denoting  BSC  noise  perturbs  this  sequence 
prior  to  being  received  at  the  decoder.   The  error  sequence 
is   (01  00  00  00  01  00  00  00  00  10  00  00  00  00  00  10  00  00 
...).   Then,  the  received  sequence  is  (10  01  00  10  11  11  00 
11  10  00  00  01  11  00  11  11  11  00  ...). 

The  following  subsections  are  the  steps  of  the  Viterbi 
algorithm  decoding  the  perturbed  (received)  sequence, 
a.   Step  1 

The  following  figure  and  tables  are  set  up  to 
indicate  the  trellis  structure,  the  path  (survivor)  sequences 
and  the  path  scores. 

(1)  Trellis.  (Figure  10). 

(2)  SSEQ(I) .   The  present  state  survivor 


sequence. 

(3)  SCORE (I) 
sequence  score. 

(4)  SSEQ(J)  . 

(5)  SCORE (J) 


The  present  state  survivor 

The  next  state  survivor  sequence. 
The  next  state  survivor  sequence 


score. 


b.   Step  2 

An  initialization  of  tables  (2)  and  (3)  is 

shown  in  figure  13.   Tables  (4)  and  (5)  are  not  considered 

until  the  first  group  of  2 (n)  received  bits  enter  the 
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decoder.   The  all  zero  sequences  placed  in  SSEQ(I)  depicts 
the  assumption  that  the  encoder  had  no  set  bits  entered 
prior  to  the  beginning  of  the  message  under  consideration. 
A  decoder  parameter  is  the  length  of  the  survivor  sequences 
used  to  hold  path  estimates  until  a  bit  is  decoded  with  the 
desired  probability  of  error.   This  length  is  called  the 
decoder  constraint  length  (DCL) ,  whose  value  is  6(2  x  K)  for 
this  example.   Also  initialized  is  SCORE (I)  so  that  a  score 
of  zero  is  in  state   0   to  denote  that  the  encoder  is 
assumed  to  have  started  from  this  state.   The  other  states 
are  assigned  scores  which  demonstrate  the  unlikelihood  of 
the  encoder  starting  from  these.   Normally,  a  value  of 
2  x  K   is  sufficient,  which  is   6   in  this  example. 
c.   Step  3  (Figure  14) 

This  step  is  the  first  in  the  actual  decoding 
process.   However,  it  will  become  obvious  that  it  is  recur- 
sive throughout  the  remainder  of  the  algorithm.   Now,  the 
first  2  received  bits  (10)  enter  the  decoder.   These  bits 
are  compared  with  each  pair  of  code  bits  on  the  branches 
of  the  trellis  and  a  distance   (A)   is  determined  for  each 
branch.   This   A   is  then  added  to  the  present  score  of  the 
corresponding  present  state.   Now  there  are  two  scores 
coming  in  on  the  branches  to  each  next  state.   At  this  time, 
a  decision  is  made  as  to  which  branch,  coming  into  each  next 
state,  has  the  lowest  score.   If  the  scores  are  identical 
then  an  arbitrary  choice  can  be  made,  decode  estimate  is 
zero  for  this  example  or  an  alternative  may  be  exercised  to 
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State 

SSEQ(J) 

Score  (J) 

0 

0 

0 
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0 
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0 
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0 

0 

0 

0 

0 

0 

6 

Decoder 
Input 

10 

Present  State  Decoded  Bit  =  0  From  State  0 


Figure  14.   Next  State  Tables  at  the  Completion  of  the 
First  Shift  of  2-Bits  Into  the  Decoder. 
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further  evaluate  the   A   added  to  the  score  in  SCORE (I). 
When  all  the  next  states  have  a  score,  they  are  entered  in 
SCORE (J).  Then  the  next  state's  survivor  sequences  are 
determined  for  entry  into  SSEQ(J) .   This  is  accomplished  by 
determining  what  bit,  0  or  1  was  dropped  from  the  present 
state  at  the  end  of  the  minimum  score  branch  for  each  next 
state.   Then,  this  bit  is  shifted  into  the  beginning  of  the 
appropriate  present  state  sequence  (SSEQ(I))  and  the  entire 
sequence  is  then  placed  in  the  SSEQ(J)  table,  corresponding 
to  its  next  state. 

The  oldest  bit  just  shifted  out  of  the  sequence  in 
SSEQ(I)  at  state   0   is  the  most  likely  decoded  bit  since 
it  is  associated  with  the  minimum  score  (0)  sequence  in  the 
present  state.   The  bit  decoding,  after  a  length  of  DCL, 
minimizes  the  effect  of  randomly  spread  errors  in  the 
received  sequence.   Again,  if  the  scores  of  two  or  more 
sequences  are  equal,  then  a  decision  is  made  according  to 
a  more  detailed  comparison  of  the  scores  in  SCORE (I). 
d.   Subsequent  Steps 

Since  the  most  recent  paths  and  scores  are  now 
in  SSEQ(J)  and  SCORE (J),  the  roles  of  present  and  next  are 
swapped,  (I)<— MJ),  and  step  3  is  carried  out  again.   Thus, 
after  each  decoding  bit  decision,  the  (I)  and  (J)  are 
swapped,  and  step  3  is  repeated  time  and  again.   Figure  15 
shows  the  decoding  steps  for  four  more  shifts  of  the  decoder 
input.   The  tables,  after  the  eighteenth  shift,  are  shown 
in  figure  16 . 
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It  is  noticed,  at  this  time,  that  the  first   8 
(DCL  +  K-l)  bits  are  zero.   This  is  explained  by  the  pre- 
sense  of  6  (DCL)  zeroes  in  the  initial  SSEQ(I)  and  the 
2 (K-l)  zeroes  in  the  encoder  when  the  first  nonzero  code 
bits  were  generated. 

Another  note  should  be  made  of  the  ease  with  which  a 
soft  decision  demodulation  scheme  could  be  incorporated 
into  the  scoring  process.   This  dimension  would  produce  a 
more  accurate  representation  of  the  likelihood  of  a  sequence 
and  possibly  eliminate  the  need  for  arbitrary  decisions, 
which  may  also  be  time-consuming. 

At  a  first  glance  this  procedure  may  seem  strange  or 
awkward,  but  there  can  be  no  doubt  that  the  Viterbi  algo- 
rithm is  easily  computer  (hardware  or  software)  implement- 
able.   The  recursive  nature  of  the  primary  decoding  step 
is  the  prime  factor  in  controlling  the  size  of  that  imple- 
mentation. 


51 


III.   COMPUTER  PROGRAM 

The  term  simulation  in  the  title  of  this  thesis  should 
not  be  construed  to  mean  the  encoder  and  decoder  operations 
(Programs)  are  simulated.   By  no  means  is  this  the  case, 
for  the  rate   1/n   convolutional  encoder  and  the  Viterbi 
decoder  are  software  implementations  that  could  be  used 
in  an  actual  system.   Of  course,  some  form  of  synchronization 
is  needed  for  practical  operation,  but  in  this  program  that 
is  assumed  to  have  been  accomplished  by  elements  (hardware 
or  software)  preceding  the  decoder  in  the  receiver. 
A.   CHANNEL  NOISE 

Actually,  the  simulation  occurs  when  a  channel  is  des- 
cribed in  the  program  by  a  noise  generating  section.   This 
section  includes  a  quantization  segment,  random  number  seg- 
ment, noise  generation  segment,  and  perturbation  (summing) 
segment.   The  following  paragraphs  are  discussions  of  these 
four  divisions  of  the  channel  noise  program. 

1 .  Quantization 

This  part  of  the  program  uses  the  encoded  sequence 
as  an  input,  then  passes  this  sequence  on  to  the  noise  sec- 
tions after   Q   zeroes  or  ones  have  been  substituted  for 
each   0   or   1.   When   Q=l   is  a  program  input  this  segment 
is  bypassed  for  apparent  reasons. 

2 .  Random  Numbers 

The  random  numbers  generated  in  this  program  were 
obtained  by  using  the  Lehmer  congruential  method. 

52 


X  .  -  =  aX   +  b  (mod  Tn) , 
n+1     n  0 

where   a  =  257,  b  =  1,  and  TQ  =  2   .   The  term  X_,   starting 

number,  is  varied  to  provide  representative  sequences  of  the 

distribution  shown  in  figure  18.   These  numbers  are  used  to 

determine  the  time  between  set  bits  in  noise  sequence.   The 

T0 
period  of  the  random  number  sequence  is   2   ,  therefore, 

T0 
a  large  sampling  of  the  sequence  approaching   2    would 

result  in  a  binomial  distribution.   However,  the  test  runs 

used  to  determine  results  in  this  thesis  employed  message 

inputs  of  10,000  bits  in  length.   Thus,  the  distribution 

(dots)  shown  in  figure  17  is  a  better  representation  of  the 

time  between  set  bits  in  the  noise  sequence. 

3 .   Noise  Generation 

After  a  random  sequence  has  been  generated  in  the 

previous  segment,   2.  ,   a  noise  sequence  is  formed  using 

the  numbers.   Multiplication  of  the  number  of  set  bits  in 

16  bits  (T) ,  the  length  of  one  random  number,  with  powers 

of   2(2   ,2   ,2  ,2  )   determine  how  many  zeroes  will  occur 

before  the  next  set  bit  in  the  noise  sequence.   For  example, 

if  the  random  number  being  considered  by  the  program  is 

0110011100100100,  then  its  weight  is  1Q   or  111,*   If  tne 

noise  parameters  (core  locations  10230  and  10236)  specify 

multiply  by  2   ,  then  the  number  7_  is  shifted  one  place 

right  and  becomes  3ft.   Likewise,  if  multiply  by  2   is 

entered,  then  70  becomes  160.   These  new  numbers  are  then 

o  o 

used  to  determine  the  space  between  set  bits  in  the  noise 
sequence.   Each  16-bit  number  from  the  congruential  random 
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number  generator  is  used  to  generate  a  space.   The  number 
of  16-bit  numbers  used  in  one  run  is  determined  by  the  num- 
ber of  code  bits  which  are  present  when  the  noise  is  added 
to  the  code  sequence . 

The  variation  in  the  set  bit  density  for  the  noise 
sequence  is  analogous  to  varying  the  density  of  channel 
errors.   This  fact  is  used  in  the  evaluation  of  the  codes 
in  section  IV.   Figures  18  through  23  show  the  probability 
of   k   errors  in  16  channel  bits  (Q  =  1,2,3)  versus  the 
number  of  errors,   k,   in  16  channel  bits.   The  plots  are 
for  averages   (XT)   of  about  10  and  2  errors  per  16  channel 
bits.   Also  shown  is  the  corresponding  plot  of  the  Poisson 
distribution  as  a  continuous  curve.   Any  discrepancy  between 
the  theoretical  and  measured  (dots)  is  attributed  to  the  16 
bit  length  constriction  primarily. 

4 .   Perturbation 

The  purpose  of  this  last  segment  of  the  noise 
simulation  is  to  simply  exclusive  OR  (additive  noise)  the 
noise  sequence  to  the  quantized  encoded  sequence.   This 
operation  forms  a  perturbed  sequence,  which  is  the  received 
sequence  for  the  decoder.   If  this  was  to  be  changed  to 
multiplication  or  some  other  operation  another  noise  form 
could  be  simulated  for  the  coding  scheme. 

All  of  the  details,  as  to  programming  specifics,  are 
discussed  in  Appendix  B,  along  with  a  program  listing. 
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Figure  19.   Simulated  Noise  (Q  =  1,  AT  =  2.0) 

vs  Poisson  Distribution   (XT  =  2.0) 
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Figure  21.   Simulated  Noise  (Q  =  2,  XT  =  2.0)  vs 
Poisson  Distribution   (XT  =  2.0) 
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Figure  23.   Simulated  Noise  (Q  =  3,  XT  =  2.1)  vs 
Poisson  Distribution  (XT  =  2.1). 
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B.   PROGRAM  FLOW 

The  four  major  sections  of  the  program  are  the  encoder, 
decoder,  noise  generator,  and  the  message  generator 
(ASCII  7-bit  code) .   An  appendix  is  devoted  to  a  detailed 
listing  and  brief  discussion  of  each  of  the  above  sections. 

The  parameters  that  are  needed  to  initiate  the  simula- 
tion are: 

1.  N;  the  inverse  rate  (1/N   ) (memory  location  10200 
in  Appendix  F) . 

2.  K;  the  encoder  constraint  length  (memory  location 
10206  in  Appendix  F0. 

3.  Q;  the  quantization  levels  (memory  location  10214  in 
Appendix  F) . 

4.  DCL;  the  decoder  constraint  length  (memory  location 
10222  in  Appendix  F) . 

5.  Noise  parameters;  operation  (-j/x)  and  operand 

(1,2 , 3 ,...)  (memory  locations  10236  and  10230 
respectively  in  Appendix  F) . 

6.  Generator  sequences;  representations  entered  so 
that  first  cell  in  encoder  is  bit   0  (memory 
locations  13760-13776,  where  G1(13760,  G2 (13762, ... , 
Gg(13776)) . 

Another  programming  concept,  position  independent  code 
(PIC) ,  is  used  to  allow  the  user  to  move  the  entire  program 
or  any  part  (with  few  changes)  to  another  part  of  core. 
This  approach  introduced  the  need  for  a  program  segment  that 
initializes  all  addresses  and  counters  (constants)  used  in 
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the  program.   The  listing  and  a  brief  description  of  this 
segment  are  included  in  Appendix  D. 

The  entire  program,  in  block  form,  is  depicted  in 
figure  24.   If  the  reader  is  interested  in  a  detailed 
program  description  he  is  referred  to  Appendices  B-E. 
Appendix  A  contains  a  map  of  that  part  of  core  used  for  the 
implementation  of  the  position  independent  code. 
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(10000-10170) 


Message  Source 
ASC  II  7  Bit  Code 


(10174-10706) 


I 


Address  and  Counter 
Initialization 


(10710-11176)^/ 


Rate  1/n  Convolutional 
Encoder 


(11200-11706 


2i 


r 


Noise  Simulator 


(11710-13176) 


\y 


Viterbi  Decoder 


(137 00-13400)  J, 


Message  Sink 
ASC  II  7  Bit  Code 


Note:   The  numbers  in  parentheses  are 
actual  addresses  occupied  by 
the  instructions  corresponding 
to  a  particular  flow  chart 
block. 


Figure  24.   Block  Diagram  of  Program  Flow. 
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IV.   RESULTS  AND  CONCLUSIONS 

The  term  best  code  is  used  in  many  papers,  but  for  the 
most  part  the  use  of  the  term  best  code  is  not  defined  in 
sufficient  detail.   This  is  one  descrepancy  that  will  be 
avoided  by  this  thesis.   Included  in  this  section  is  a 
summary  of  the  computer  results  obtained  from  a  best  code 
(defined  below)  determination  procedure  and  encoder/decoder 
parameter  variations. 

A.   BEST  CODE  DETERMINATION 

The  determination  of  a  best  code  for  a  specific  channel 

a  simulated  channel  in  this  paper,  follows  a  procedure  based 

on  Shannon's  fundamental  theorem  for  a  discrete  channel  with 

noise.   This  procedure  is  best  described  using  Shannon's 

representation  of  the  attainable  region  in  a  graph  (figure 

25)  of   H(x)   (information  rate)  versus   H  (x)   (probability 

of  decoded  messager  error).   If   H(x)  <  C   (channel  capacity), 

then  Shannon  shows  that   H  (x)   can  be  made  arbitrarily  small 

Y 

with  a  proper  encoding  procedure.   When   H(x)  >  C,   then  the 
excess  information  being  pushed  onto  the  channel  can  only 
increase  the  uncertainty   (H  (x) )   of  the  decoded  message. 
The  minimum  value  of   H  (x)   is  very  close  to   H-C   in  this 
latter  case. 

The  computer  program  in  this  thesis  is  implemented  with 
a  supplemental  segment  which  stepped  through  all  well-defined 
convolutional  code  generators  of  a  specific  code  rate  and 
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constraint  length.   Well  defined  refers  to  the  constraint 
length  definition  of  the  code,  meaning  that  the  coefficient 
of  the  zero  power  in  the  generator  polynomial  must  be  set  for 
at  least  one  generator  sequence,  and  then  another  sequence 
must  have  the  coefficient  set  for  the   K-l   power.   An 
example  of  this  is  shown  here: 

n  =  2,  K  =  4;   2X (D)  =  1  +  D2 

g_2(D)  =  D  +  D2  +  D3 

n  =  3,  K  =  3;   g,  (D)  =  D  +  D3  =»  1  +  D2 

g_2  (D)  =  D2  +  D3  =»  D  +  D2 

g_3(D)  =  D  +  D2  +  D3=^  1  +  D  +  D2 

The  rate  is  defined  by  the  number  of  nonzero  generator 
sequences. 

The  first  two  paragraphs  describe  a  procedure  and  the 
selection  of  encoders  to  be  processed  by  the  program.   The 
program  is  applied  to  find  the  generator  sequences,  for  a 
given  code  rate/constraint  length,  which  come  the  closest  to 
the  lower  error   (H  (x) )   boundary  of  figure  25.   Code  rates 
of  1/2  and  1/3  are  tested  and  if  the  input  message  is 
assumed  to  have  1  bit  of  information  in  each  bit  entering 
the  encoder,  then  these  rates  correspond  to  those  plotted 
on  the  horizontal  scale  of  figure  25.   Tables  II  through  IV 
present  a  summary  of  the  test  run  results  for  a  best  code 
determination. 

Each  test  run  was  made  for  a  40,000  bit  encoder  input 
applied  over  various  samples  of  the  simulated  noise  and  at 
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four  values  of  XT(1,2,4,10  errors  per  16  bits).  The  com- 
puter output  is  a  listing  of  the  decoded  message  errors  for 
each  well  defined  encoder.  This  output  would  require  many 
additional  pages  to  list  in  this  thesis.  Therefore,  the 
summary  in  Tables  II  through  IV  is  presented  instead.  The 
number  of  decoded  errors  is  then  divided  by  4  0,000  to  give 
the  measure  of  average  error  for  each  test  run   (XT) . 

The  various  values  of   XT   are  obtained  by  the  noise 
simulation  and  have  the  characteristic  distributions  of  those 
curves  presented  in  figures  18  and  19. 

A  discussion  of  these  results  is  presented  in  part  C  of 
this  section. 

B.  DECODER  PARAMETER  VARIATIONS 

In  order  to  provide  an  insight  as  to  the  effects  of 
quantization  (Q)  and  decoder  constraint  length  (DCL) ,  tables 
V  and  VI  are  shown.   The  effect  of  increasing  Q  or  DCL  is  to 
decrease  the  frequency  of  errors  thus  improving  the  error 
correcting  scheme.   The  codes  used  in  these  tables  are  those 
chosen  as  the  best  codes  from  tables  II  through  IV. 

C.  DISCUSSION  AND  CONCLUSIONS 

The  determination  of  a  best  coding  scheme  for  a  specific 
channel  application  can  be  made  from  tables  II-IV.   The 
generator  sequences  which  have  the  minimum  frequency  of 
errors  (in  place  of   H  (x) )  for  the  various  values  of   XT 
are  outlined  boldly.   Using  this  comparison  technique  plus 
a  comparison  of  the  general  statistics  of  the  encoder  at 
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other  values  of  XT,  the  following  generator  sequences  are 

chosen  to  be  the  encoders  which  meet  the  criteria  set  forth 
previously. 

N  =  2,  K  =  3;   2l  =  (110),  g_2  =  (001) 

N  =  2,  K  =  4;   2l  =  (1110),  g^  =  (0101) 

N  =  3,  K  =  3;   2l  =  (111),  g_2  =  (001), 

£3  =  (001) 

These  encoders  were  also  used  for  the  results  obtained  by 
varying  Q  and  DCL  in  tables  V  and  VI. 

With  actual  concrete  measurements  of  the  effectiveness 
(minimum  error)  of  a  code  and  a  means  (computer  program)  of 
applying  these  codes  to  actual  channel  recordings  there  is 
no  reason  to  have  to  settle  for  an  inadequate  error  correcting 
system.   The  results  presented  in  this  section  indicate  clear 
preferences  in  choosing  certain  codes  to  accomplish  desired 
communications  in  a  given  noise  environment.   The  user  speci- 
fics the  task  and  with  the  computer  a  code  can  be  chosen. 
The  true  value  of  such  a  program  should  be  evaluated  when 
implemented  with  actual  channel  noise.   No  one  has  been  able 
to  find  a  relation  or  algorithm  for  encoding  that  would  enable 
a  communication  system  to  reach  the  maximum  average  informa- 
tion rate  (c)  with  an  arbitrarily  small  frequency  of  error. 
However,  computer  analysis  offers  the  channel  user  an  oppor- 
tunity to  find  a  coding  scheme  which  meets  the  requirements 
for  that  channel. 
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APPENDIX  A 
POSITION  INDEPENDENT  CODE  CORE  MRP 

THE  CORE  LOCRTIONS  OF  ADDRESSES,  COUNTERS,  RND 
MISCELLRNEOUS  STORRGE  (WORK)  AREAS  ARE  LISTED  BELOW. 
SOME  OF  THESE  VALUES,  PRECEDED  BV  *,  ARE  ENTERED  BV 
OR  COMPUTED  BV  THE  PROGRAM  INITIALIZATION  SUBPROGRAM, 
WHICH  FOLLOWS  THE  INPUT  SUBPROGRAM. 

*  15660;  INPUT  BLOCK  RDDRESS 

*  12682;  ENCODED/DECODED  BLOCK  RDDRESS 

*  12604;  QUANTIZED  CODE  SEQUENCE  BLOCK  ADDRESS 

*  12606;  NOISE  SEQUENCE  BLOCK  ADDRESS 

*  12610;  RANDOM  NUMBER/PERTURBED  SEQUENCE  BLOCK  RDDRESS 

*  12612;  TRELLIS  TABLE  BLOCK  ADDRESS 

*  12614;  S5EGKI)  TABLE  BLOCK  RDDRESS 

*  12616;  SCOREU)  TABLE  BLOCK  ADDRESS 

*  12620;  SSEQCJ)  TABLE  BLOCK  ADDRESS 

*  12622;  SCORE <J>  TABLE  BLOCK  ADDRESS 
12624;  "NOT  USED" 

*  12626;  GENERATOR  SEQUENCE  BLOCK  ADDRESS 

*  12620;  STACK  ADDRESS 
12622;  HOLD  (STORRGE) 

12  624;  D  E  L  T  A  <  U  S  E"  D  I  N  D  E  C  0  D  t  R  ) 

12626;  MfiSK  (USED  TO  MAKE  DECODED  BIT  DECISION) 

12640;  WORK  (STORAGE) 

12642;  MESSAGE  WORD  <16  BITS)  COUNT  WORK  LOCATION 

12644;  MQ  BIT  COUNT  WORK  LOCATION 
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13646; 
13650; 
13652; 
13654, 
13656; 
13660, 
1366; 
13664, 
13666, 
13670; 
136?; 
13674, 
13676; 
+   13700, 

*  i37o; 

*  13704, 

*  13786 
13718; 
13712; 
13714. 
13716, 
13728, 
1372, 
13724; 
13726, 
13730, 
13732; 


DECODED  BIT  COUNTER 

NUMBER  OF  SCORES  COUNTER 

MINIMUM  SCORES  ADDRESS 

MINIMUM  SCORE  VALUE 

N-BIT  LOOP  COUNTER 

ZERO  BIT  SCORE  VALUE 

ONE  BIT  SCORE  VALUE 

SSEQ  WORK  LOCATION 

TOTAL  OF  NOISE  VALUES  IN  NOISE  SEQUENCE 

TOTAL  NUMBER  OF  ERRORS  IN  NOISE  SEQUENCE 

TOTAL  NUMBER  OF  DECODED  MESSAGE  BIT  ERRORS 

"NOT  USED" 

"NOT  USED" 

NUMBER  OF  CODE  BITS  PER  INPUT  BIT,  N 

ENCODER  CONSTRAINT  LENGTH,  K 

QUANTIZATION  VALUE,  Q 

DECODER  CONSTRAINT  LENGTH,  DCL 

MESSAGE  WORD  C16  BITS)  COUNT,  W 

POWER  OF  TWO  USED  IN  GENERATING  NOISE  SEQUENCE 

0  P  E  P  f  I T I  0 N  <  D I  V I D  E / M  U L  T I P  L  V )  USE  D  F  0  R  NOISE 

Q  X  N 

Q  X  N  X  W 

4  X  K  (NONZERO  STATES  INITIAL  SCORE) 

» 

2  **  K 

<  2  *  *  K )  -  1 
2  **  <K  -  1) 
<2  **  <K  -  1))  -  1 
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* 

12724; 

* 

12726; 

* 

12740; 

* 

12742; 

* 

12744; 

* 

12746; 

* 

12756; 

12752; 

12754; 

12756; 

12760; 

12762; 

12764; 

12766; 

12770; 

12772; 

12774; 

12776; 

2  **  <K  -  2) 

<2  **  G!)  -  1 

DCL  +  K  -  1 

2  *.*  <DCL  -  1) 

SCORE<I>  -  5SEGKI)  CSUBTRRCT  ADDRESSES) 

SSEGKJ)  -  SSECKI)  (SUBTRRCT  ADDRESSES) 

n  x  m 

NEXT  SCORE  ADDRESS 

NEXT  SSEQ  ADDRESS 

STORAGE  LOCATION  FOR  A  COUNTER  VALUE  OF 

G<1)<  GENERATOR  SEQUENCE 

G<2),  GENERATOR  SEQUENCE 

G<3)j  GENERATOR  SEQUENCE 

G<4)i  GENERATOR  SEQUENCE 

G<5),  GENERATOR  SEQUENCE 

G<6>,  GENERATOR  SEQUENCE 

G<?>i  GENERATOR  SEQUENCE 

G<8)i  GENERATOR  SEQUENCE 


THE  WORD  COUNT  IN  LOCATION  12719  IS  PLACED  THERE  BV 
AN  INSTRUCTION  AT  THE  END  OF  THE  MESSAGE  INPUT  SUB- 
P R 0 G R Fl M .   T H E  V A L U E  OF  "i  H E  G E N E R  A T  0 R  S E 9. U E N C ES  I  S 
E  N 1  E  R  E  D  B  V  7  H  E  P  R  0  G  R  Fl  M  N  E  R  B  E  F  0  R  E  S  T  Fl  R  T  I  N  G  A  R  U  N . 
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APPENDIX  B 
PROGRAM  LISTING  OF  NOISE  SIMULATION 

The  following  flow  chart  depicts  the  flow  of  the 
instructions  in  the  following  machine  language  listing 
of  the  noise  simulation  subprogram.   The  numbers  (base  8) 
on  the  upper  left  of  the  blocks  in  the  flow  chart  corres- 
pond to  those  instruction  addresses  of  the  subprogram 
mentioned  in  that  block. 
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(11200 

-11324) 

This 

segment 

quantizes  the 

encoded 

block 

,  that  is,  encoded  1' 

s  become 

Q  l's 

and  encoded  0 ' s  become  Q  0 ' s . 

The  result  is  then  placed 

in  block 

speci 

fried  at 

(13604)  . 

(11330-11374) 


1 


This  segment  generates  the  random 
number  sequence  by  the  Lehmer  Con- 
gruential  Relation.   These  numbers 
are  stored  in  a  block  beginning  at 
the  address  specified  at  (13610) . 


(11376-11546) 

This  segmen 

t  uses  the 

noise  para- 

meters  plac 

ed  in  locat 

10ns 

,  (10230) 

and  (10236) 

to  change 

the 

value 

of  the  number  of  set  £ 

>its 

in  a 

16 

bit  random 

number,  so 

that  the 

density  of 

errors  (AT) 

in 

the  noise 

is  varied. 

The  noise 

sequence 

is 

stored  in  a 

block  beginning  at 

an 

|  address  in 

(13606) . 

(11550-11646)      I 
This  segment  determines  the  noise 
value  total  (13666)  and  the  number 
of  errors  (13670)  for  a  given 
message,  i.e.,  a  value  ^  4  is  an 
error  for  Q=3. 


(11650-11700) 


I 


This  segment  adds  the  noise  block 
to  the  quantized  block  and  stores 
the  result  in  a  perturbed  block 
beginning  at  an  address  in  (13610) 

>' 

END  OF  SUBPROGRAM 
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The  following  computer  printout  is  the  noise  program 
just  discussed. 


011268 
611202 
011204 
011206 

011210 

011212 
011214 
011216 
011220 
011222 
011224 
011226 
011230 
011232 
011234 
011236 
011240 
01124  2 
01124  4 
011246 
011250 
011252 
011254 
011256 
011260 
011262 
011264 

01126  6 
011270 
011272 
011274 

01127  6 

011300 

0113 0  2 
011364 
0113  0  6 
011 3 1 0 
011312 
011314 
011316 
011320 
011322 
011324 
011326 


•016760 
• 0  0  2  3  7  6 
/  0 1 6  7  0 1 
• 0  0  2  3  7  4 
•016702 
/  0  0  2  5  3  4 

•  0  2  2  7  6  7 

• 0  0  0  0  01 

•002462 

• 0  010  0  3 

•  6 1  2  0  2 1 
• 0  7  7  2  0  2 

•  0  0  0  4  3  6 
• 012  70  3 

•  0  0  0  0  2  0 

•012704 

•  0  0  0  0  2  0 

•016705 

•  0  0  2  4  3  6 

•  0  0  6  3 1  0 

•  1  0  3  4  6  5 

•  0  0  6  311 
•005303 

• 0  01410 

•  0  7  7  5  0  4 

•  0  0  0  41 6 

•  0  0  6  311 
•005211 

•  0  0  5  30  3 

•  0  014  0  6 
• 0  7  7  5  0  5 

•  0  0  0  410 

•  0  6  5  7  21 

•  012  7  0  3 

•  0  6  8  0  2  0 
•000764 

•  0  0  5  7  21 

•  012  7  0  3 

•  0  0  0  0  2  0 

•  0  0  0  7  6  6 

•  0  7  7  4  3  0 

•  0  0  5  7  2  0 

•  0  7  7  2  3  4 

•  0  0  0  0  0  0 


011330 
011332 
011334 
611336 
011340 
01134  2 
011344 
611346 
011350 
011352 
011354 
011356 
011360 
611362 
611364 
01136  6 
011370 
011372 
011374 
011376 

011400 

011402 
0114  0  4 


•01670 0 

•002254 
•012701 
•004704 

•012702 

•  0  0  0  4  01 

•012737 
•044444 
•177304 
•010237 
•177306 
•605237 
•177304 

•  0 1 3  7  2  0 
•177304 
•077107 
•010067 

•0  0230  0 
•00  024  0 

•005067 

•  0  0  2  2  6  4 

•  0  0  5  0  6  7 

•  0  0  2  2  6  2 


011406  /'016700 
011410  •002176 


011412 
011414 
011416 
011420 
011422 


• 016  7  01 

•  0  0  217  0 

•  012  7  0  3 

•  0  0  0  0  2  0 

•012704 


0114  24  ,-'000020 
0114  26  ,-'012006 


6114  3  0 
0114  32 
011434 
011436 
0114  4  0 
6114  4  2 
011444 
011446 
0114  5  0 
011452 
0114  5  4 
0114  5  6 


•020067 
•002240 

•  0  0 1 0  0  3 

•  0 1 6  7  0  0 
• 0  0  214  6 
/  0  1  2  0  0  6 

•0  0500  5 

•  0  0  6  2  0  6 
• 0  0  5  5  0  5 
• 0  7  7  4  0  3 

•  0  0  5  7  6  7 

•  0  0  2  2  3  4 
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011466 
611462 

611464 
611466 
011470 
011472 
011474 
011476 

011560 

011502 
011504 
01150  6 

01151O 

011512 
011514 
011516 
011520 

r4  4  A  C  "x  o 
Oil  J&.C 

0115  2  4 
011526 
011538 
011532 

011534 
011536 
011540 
01154  2 
811544 
011546 
8115  5  0 
811552 
811554 
811556 
8 115  6  0 
811562 
8115  6  4 
8115  6  6 
8115  7  8 
811572 


/ 8  0 10  8  5 

/816704 
/  0  0  2  2  2  4 
700  6  20  5 
7  0  7  7  4  0  2 

7 0  0  0  4  0  7 

7016737 
7  0  0  2  2 1 2 
7177304 

/  010537 

7 17  7  2  0  6 
7 013  7  0  5 
,-'17  730  4 
/  0  0  5  3  8  3 
7  0  0  2  4  0  5 
7  0  0  6  3  1 1 
/  0  0  5  S  0  5 
7  6  B  3  3  ?  3 
/  8  0  5  211 
7  0  0  0  7  2  5 
7 0  05 7  21 
7 0  2 016  7 
7  0  0  2  0  5  2 

7  0  0  10  01 
/ 0  0  0  4  0  2 

7012703 

7 0  0  0  017 

7  0  0  0  7  6  2 
7  016  7  0  0 
7  0  0  2  0  3  2 
7  0 1 6  7  01 
,-'002140 
/ 012  7  0  2 

7  0  0  0  0  2  0 

7  0 1 2  0  0  4 
7  8 1 6  7  0  3 
7 0  0  211 2 
/  0  0  5  0  0  5 


811574 
011576 

011600 

01160  2 
011604 
01160  6 
811610 
811612 
011614 
011616 
811620 
811622 
011624 
811626 
8116  2  8 
811632 
811624 
811626 
811648 
011642 
811644 
811646 
811658 
811652 
811654 
811656 
811668 
811662 
811664 
8116  6  6 
811678 
811672 
811674 

8116  7  6 
811780 

8117  8  2 
* 


Z006205 
7806204 
/ 0  0  5  5  0  5 
7005202 
7  0  0141 2 
,•'077306 
7060567 
7  8  0  2  0  5  2 
7  016  7  0  2 
7  0  0  2116 
7  0  0  6  2  0  2 
7  0  2  0  5  0  2 
7  0  0  2  7  6  0 
7005267 
7002826 
7000755 
7  012  0  0  4 
7  812  7  0  2 

70  0  0  820 

7077117 
7  0  0  0  2  4  0 
7  0  0  0  2  4  0 
7016700 
70  01732 
7016781 
7  0  017  2  4 
7  016  7  0  2 
7  0  017  2  4 
7  0 1 6  7  0  5 
7  0  0  2  0  3  0 
7  012  8  0  3 
7  01210  4 
7  0  7  4  2  0  4 
7  010  4  2  2 
,-'  0  7  7  5  0  5 
7  £l  0  0  0  0  0 
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APPENDIX  C 

Program  Listing  of  Rate  1/n 
Convolutional  Encoder 

The  following  flow  chart  depicts  the  flow  of  the 
rate  1/n  convolutional  encoder  implementation.   The  cor- 
responding machine  language  program  follows  the  flow 
chart  and  the  instruction  addressed  are  indicated  on  top 
of  the  respective  flow  chart  block  to  which  they  corres- 
pond in  the  subprogram. 


-> 


(10710-10726) 


This    £ 
of   the 
coded 
(R6) 

>egment    sets   up 
:   message   block 
block    (Rl) ,    and 

the 
(R0) 
the 

addresses 
,    the   En- 
stack 

(10730- 

•10772) 

•V 

This  segment  uses  the  extended 
arithmetic  element  (AC-MQ)  to 
shift  in  a  new  message  bit  to  be 
encoded  with  the  previous  entries 
within  the  encoder  constant  length, 


(10774-11062) 


vy 


This  segment  determines  the  bits  in 
AC (177302)  that  are  to  be  encoded 
and  steps  through  the  generator 
sequence  block  computing  the  n  out- 
put  code  bits. 


(11064-11110) 


^ 


This  segment  contains  loop  indica- 
tors for  the  number  of  bits  left  in 
MQ  and  the  message  word  count.   If 
the  latter  is  not  exhausted,  a  loop 
is  made  to  put  a  new  16-bit  word  in 
MQ (177304)  . 


END  OF  SUBPROGRAM 
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The  following  computer  printout  is  the  encoder  implementa- 
tion just  discussed. 


610710 

•016700 

610712 

•  0  0  2  6  6  4 

010714 

•  0 1 6  7  0 1 

610716 

•  0  0  2  6  6  2 

010720 

•012704 

010722 

/  0  0  0  0  2  0 

010724 

/  016  7  0  6 

010726 

•002706 

010720 

•016716 

010732 

•  0  0  2  7  5  4 

B10734 

•012037 

010736 

/177304 

010  74  0 

•005037 

010742 

•177302 

010  74  4 

•  0  0  0  4  0  6 

010  7  4  6 

•  0 1  3  7  4  6 

010750 

•  1  7  7  3  0  2 

010752 

•012037 

010754 

•177304 

010756 

•012637 

010760 

•  1  7  7  3  0  2 

010762 

•  0 1 2  7  4  6 

010764 

•  0  0  0  0  2  0 

010766 

•012737 

010770 

•  0  0  0  0  0  1 

010  7  7  2 

•  1 7  7  3 1 4 

010774 

•01674 6 

010776 

•  0  0  2  7  0  0 

011000 

•  0  1  3  7  0  5 

0110  0  2 

•177302 

0110  0  4 

•016702 

0110  0  6 

•  0  0  2  61 6 

011010 

• 0  0  510  5 

011012  •012302 

611014  ^040502 

611016  •005046 

611020  •016746 

011022  •002656 

011024  •005726 

011026  ^006202 

611030  ^005516 

011032  •685346 

011034  •061373 

611036  •ooesn 

611040  •665726 
0110  4  2  • 0  0  6  2  2  6 
611044  •005511 
011046  ,-'005304 

01105  0  •001003 

011052  •012704 
011054  •888828 
011056  •005721 
011060  ,''005316 
011062  •001353 
011064  ,-'005726 
011066  /&B5216 
011070  •001336 
011072  •665726 
611074  ,-'005316 
011076  ,-'001323 

01110  0  • 0  0  5  7  0  4 

01110  2  • 0  0 1 4  0  3 

011104  •006311 

01110  6  ,-'  0  0  5  2  0  4 

011110  .'001375 

011112   Z000000 
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APPENDIX  D 

Program  Listing  of  Viterbi  Decoder 

The  following  flow  chart  depicts  the  flow  of  the 
Viterbi  decoding  algorithm  implementation.   The  correspond- 
ing machine  language  program  follows  the  flow  chart  and  the 
instruction  addresses  are  indicated  on  top  of  the  respec- 
tive flow  chart  block  to  which  they  correspond  in  the 
subprogram. 
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(11710-12042) 


COMPUTE  TRELLIS 
AND  STORE  IN  TRELLIS 
BLOCK  (13612) . 


(12044- 

-12146) 

i 

QUANTIZE  TRELLIS 
ENTRIES  FOR   Q 

TABLE 
>  1 

(12150- 

-12206) 

7 

INITIALIZE  SSEQ(I)  (13614)  AND 
SCORE  (I)  (13616) 


(12210-12366)^ r 


DECODE  OLDEST  BIT  IN  THE 
SSEQ(I)  SEQUENCE  WITH 
THE  MINIMUM  SCORE 


(12370-12546) 


I 


SHIFT  INTO  MQ (177304) 
N  BITS  FROM  THE  PERTURBED 
SEQUENCE  (13610) 


(12550-12772) 


COMPUTE  A  BETWEEN  THE  N-BIT  INPUT 

AND  EACH  N-BIT  BRANCH  CODE.   ADD 

THE  BRANCH  A  TO  IT  CORRESPONDING 

PRESENT  SCORE.  DETERMINE  THE 

BRANCH  WITH  THE  MINIMUM  SCORE 

ENTERING  EACH  NEXT  STATE. 


(12774-13072)    i 

CHANGE  EACH  SSEQ 

CORRESPONDING  TO  THE 

LOST  FROM  THE  PRESENT 

BIT 
STATE 

(13074-13166) 
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The  following  computer  printout  is  the  decoder  implementa- 
tion just  discussed. 


en  ri& 

011712 
611714 
011716 
011720 
011722 
011724 
011726 
01172  0 
011732 
011724 
011726 
011740 
011742 
011744 
011746 
011750 
011752 
011754 
011756 
011766 
011762 
011764 
011766 
011770 
011772 
GU1774 
011776 

012000 

0120  0  2 
01200  4 
012006 

012010 

012012 
012014 
GU2016 
012028 
012022 
012O24 
01202b 
0120  20 
012022 
0120  24 
0120  2  6 
012040 
012042 
0120  4  4 
0 12  0  4  6 


/016?00 
•002014 
•816701 

•001706 
•016702 
• 0  016  6  6 
• O16  70  4 
• 0  01?  0  0 
/0 0500 2 

•005012 
•016714 
•00174 0 

• 0  0  51 0  2 

•  0 1 2 1 0  5 

•  0  4  0  2  0  5 
•016744 
•0  01720 

•  0  0  5  0  0  6 
•0  06205 
•005506 
• O05214 

•  0  0 12  7  4 

•  0  0  6  2 1 2 
•006206 

•  0  0  5  51 2 
•005724 
• 0  0  5  214 

•  0  0 12  61 

•  0 1 6  7  0 1 

•  0  0 1 6  2  2 

•  0  0  5  7  2  2 

•  0  0  51 0  2 

•  0  2  6  7  0  2 
• 0  01712 
• 0  014  0  7 

•  0  2  0  2  0  0 

•  0  0  2  0  0  2 

•  0  6  0  0  0  2 

•  0  0O74  2 

•  1 6  0  O  0  2 

•  0  6  5  2  0  2 

•  0  0  0  7  2  7 

•  0  2  2  7  6  7 

•  0  0  0  0  0 1 

•  0  016  4  2 

•  0  0  1  4  4  2 

•  016  7  0  0 

•  0  015  4  2 


012050 
012052 
012054 
012056 
0120  6  0 
0120  6  2 
012O64 
012066 
012070 
01207  2 
012074 
012076 

012100 

012102 
012104 
012106 

012110 

012112 
012114 
012116 
012120 
012122 
012124 
012126 
012120 
012122 
012124 
012126 
012140 

01214  2 
012144 
012146 
012150 

01215  2 
01215  4 

01215  6 

01216  0 
812162 
01216  4 
012166 
012170 
012172 
012174 
012176 
012  20  0 
0 1  2  2  0  2 
0 1 2  2  0  4 
O1220  6 


•  0 1 6  7  0 1 

•  0  0 1 6  5  0 

•  016  7  0  2 

•  0  0 1 5  5  0 

•  0 1 6  7  0  2 

•  0  0 1 6 1 4 
• 016  7  0  4 

•  0  0 1 6  1  4 
•005012 

•  0  0  6  2 1 0 
•102402 

•  0  0  6  2 1 2 
•077402 

/ 0  0  0  4  0  2 

•0  06212 
•005212 
•077402 
•005742 
•077215 

•  0  0  5  7  2  2 

•  0 1 6  7  0  2 
•0  01554 

•  0  0  5  01 0 

•01670  4 

•  0  015  5  2 
•006210 

•  0  7  7  4  0  2 

•  0  6  2  2 1  0 

•  0  7  7  3  0  6 

•  0  0  5  7  2  0 

•  0  0  5  7  4  2 
•077126 
•016700 

•  0  0 1  4  2  6 

•  0 1 6  7  01 

•  0  014  2  4 

•  0 1 6  7  0  2 
001422 
O16702 
O014  2  0 
0  0  5  0  21 

•005O22 

•  0 1 6  7  0  4 
•0O1522 

•  0  0  5  0  2 1 

•  0 1 6  7  2  2 
•001514 

•  0  7  7  4  0  4 


86 


012210 
012  212 
012214 
012216 
012226 
012222 
012224 
012226 
012228 
012222 
012224 
01222  6 
012240 
012242 
012244 
012246 
012258 
012252 
012254 
012256 
012260 
012262 
012264 
01226  6 
012270 
012272 
012274 
012276 
0 1 2  2  0  0 
01220  2 
012204 
012206 
012210 
012212 
012214 
012216 
0 1 2  2  2  0 
012222 
012224 
012226 
012  2  2  0 
012222 
012224 
012236 
012340 
0 12  3  4  2 
012244 
012  3  4  6 


•016701 

•  0  0 1  4  0  0 
•016702 

/  0  0 1  3  7  6 
•016704 

•  0  0 13  5  6 

•  016  7  6  7 

•  0  0 1  4  6  0 

• 0  01410 

•  012  7  6  7 

•  0  0  0  0  2  0 

•  0  0 1 4  0  4 

•  012 3  3  7 
•177304 
•01276  7 

•  0  0  0  0  2  0 
/ 0  012  7  4 

•  0  0  5  2  6  7 
•001462 

•  0 0  2  0  51 
•016767 
• 0  014  4  4 
•001262 
•010267 

•  0  0 1 2  6  0 
•012267 

•  0  012  5  6 

•  0  0  5  2  6  7 

•  0  012  4  6 
•001405 

•  0  2 12  6  7 
• 0  012  4  4 

•  0  0  2  7  6  6 

•  005  722 

•  0  0  0  7  7  0 
•166767 

•  0  014  2 2 

•  0  0 1  2  2  6 
•016706 

•  0  0 1  3  2  2 
•006314 

•  0  3  6  71 6 
/  0  0 1  4  0  4 

•  1  0  0  0  0  4 

•  e  b  0  2  4  0 

•  0  0  0  2  4  0 

•  0  0  0  2  4  0 
• 0U5214 


012350 
012252 
012254 
012256 
012  36  0 
012362 
012264 
012266 
012270 
012272 
012274 
012276 
012400 
012402 
0124  04 
0124  0  6 
012410 
012412 
012414 
012416 
012420 
012422 
012424 
012426 
01242  0 
012422 
012424 
012426 
012  4  4  0 
012442 
0124  4  4 
0124  4  6 
012450 
012452 
012454 
012456 
O12460 
012462 
012464 

0124  6  6 
GU2470 
012472 
012474 
012  4  7  6 
0 1 2  5  0  0 
0125O2 
012504 

0125  06 


•005267 
•001272 

•001012 

•005724 
•005267 
•001256 

•  0  0 1 0  0  2 

•000000 

•  0  0  0  2  4  0 

•  0  0  0  2  4  0 
•012767 

•0  0  0020 

•001244 

•010102 

•066  702 

•  0  0 1  3  3  4 

•  0  0  5  0  3  7 

•  1 7  7  2  0  2 
•016706 
•601276 
•012727 

•000  001 

•177214 

•  0  0  5  3  6  7 

•  0  0 1  2 1  2 
•00140  2 
/&77S&? 
• 0  0  0  414 

•  012  7  6  7 

•  0  0  0  0  2  0 

•  0  0117  6 
•012767 
•177202 

•  0  0 1 1  5  6 
•012227 
•177204 
•016727 
•001146 

•  1  7  7  3  0  2 

•  0  0  0  7  6  2 

•  0 1 6  7  6  7 

•  0  0124 0 

• 0  011 6 0 

•  0 1 0  2  6  7 

•  0  0 1 2  5  0 

•  0 1  0 1  6  7 
•00124 6 

•  0  2  0  2  6  7 
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612516 

612512 

6125-14 

612516 

612520 

612522 

612524 

612526 

61253:0 

012512 

612534 

012536 

612540 

012542 

012544 

012546 

612550 

012552 

012554 

612556 

6125  6  6 

612562 

612564 

612  5  6  6 

612570 

012572 

612574 

6125  7  6 

0126  0  6 

O12602 

012604 

0126  0  6 

012616 

012612 

012614 

612616 

612620 

012  6  2  2 

612624 

0126  2  6 

012  6  2  0 

612632 

012634 

0126  2  6 

01264  0 

01264  2 

012644 

0 1 2  6  4  6 


•  0  0  1  1  6  6 
•100010 

•  0667  6  7 
•661226 
•061236 
•666767 
•661226 
•661224 
• 0  6  0  4  6  7 

•  0  0  0  2  4  0 
/ 1 6  6  7  6  7 

•  6  012  0  6 
•001210 
/ 1  6  6  7  6  7 
•001200 
/  0  0 1 2  0  4 

•  012  7  6  7 
/  0  0  0  0  6  2 

• 0  012  0  0 

•011267 

/  0  0 1 0  7  6 
/  0 1 6  7  0  6 

•001112 

•013767 

,i  a  -?  y  -■  r i  3 
<  J.  i  i  5  O  c 

•  0  010  4  4 
•012005 

•  0  7  4  5  6  7 
•001036 

•  0  0  5  0  6  7 

•  0  0 10  2  6 

•  O  0  5 1 6  7 
• 0  610  2  6 

•  0 1 6  7  0  5 

•  0  0 1  0  6  6 

•  016  7  6  7 

•  0  01114 
/ 0  0101 2 

•  0  4  6  7  6  7 

•  0  01016 

•  0  0 10  0  4 

/&€€?£■? 

•  0  O  1  0  0  0 

•  0  0  0  7  7  4 

•  0  O  6  2  6  7 

•  0  0O774 
• 0  7  7  5  0  3 
•077617 


012656 

• 0  0  0  24  0 

612652 

•666767 

612654 

•  0  0  6  7  5  6 

612656 

•  0  0  1  0  0  0 

6126  6  0 

•010267 

612662 

•  0  0  6  7  4  6 

612664 

•  0  6  6  7  6  7 

612666 

•  0  010  4  0 

612670 

•000740 

612672 

•617767 

612674 

•  0  6  0  7  3  4 

612676 

•  0  0  6  7  6  2 

61276  0 

•016706 

612762 

• 0  0  6  7  7  4 

612764 

•  0 1 3  7  6  7 

612766 

•177302 

612716 

•  0  0  O  7  2  6 

612712 

•012005 

612714 

•  0  7  4  5  6  7 

612716 

•000726 

012726 

•  0  0  5  0  6  7 

612722 

•  0  0  0  7 1  0 

612724 

•  0  0  5 1  6  7 

612  7  2  6 

•000710 

612736 

•616705 

612732 

• 0  0  0  7  5  6 

012734 

• 016  7  6  7 

012  7  3  6 

•  0  &&7  7 6 

012740 

•0O0  674 

012742 

•046767 

612744 

•000672 

612746 

•  0  6  6  6  6  6 

6127  50 

• 0  6  6  7  6  7 

012752 

•  0  0  0  6  6  2 

012754 

•  6  0  0  6  5  6 

0127  5  6 

•  0  0  6  2  6  7 

O12760 

•  0  0  0  6  5  6 

6127  6  2 

•  0  7  7  5  0  3 

61276  4 

•  0  7  7  617 

612766 

•  &  6 67 6  7 

012776 

•  O  O064  2 

612772 

•  0  O  O  6  6  6 

012774 

•  0  2  6  7  6  7 

012776 

•  0  0  0  6  6  0 

0  1 3  O  0  0 

•  0  0  0  6  6  0 

01300  2 

•  0  0  3  4  21 

0130O4 

•  016  7  0  5 

01300  6 

•  0  0  0  7  4  2 
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013EU0 

/  9 1  6  7 1  5 

612012 

Z000646 

012014 

/  0 1 8 1 6  7 

012016 

/ 090644 

012020 

/Q6676? 

01202  2 

/ 0  0 0  7  0  4 

012O24 

/ 0  0  @63  6 

012026 

/  0 1 6  7  0  5 

0120  2  0 

/  8  0  0  6  3  2 

012O22 

/  0 1 6  7  0  6 

012O34 

,'  0  0  0  716 

012O3  6 

,-'011516 

012049 

/ 0  0  6  216 

012042 

/  0  0  5  2 1 6 

012044 

/ 0  0  0  41 0 

012046 

/016705 

01205  0 

/  0  0  0  7  0  0 

0120  5  2 

/  0 1 6  7  0  6 

012054 

/  0  6  0  6  7  6 

012056 

,-■016  715 

01206  0 

/ 0  0  0  5  7  6 

012O6  2 

/011116 

012064 

/006316 

012O66 

/016705 

012O70 

/  0  0  0  6  6  0 

0 1 2  0  7  2 

,-022526 

CH2074 

/  0 1 0  5  6  7 

0120  7  6 

1 0  0  B  6  5  2 

013100 

012102 
012104 
012106 

012110 

012112 
012114 
012116 
012120 
012122 
012124 
012126 
012120 
012122 
012134 
012126 
012140 
012142 
012144 
01214  6 
012150 
012152 
012154 
012156 
012160 
012162 
012164 
01216  6 


/  0 1 0  6  6  7 

/ 0  0  0  6  5  0 

/ 0  0  5  2  6  7 
,-'0  0  0  64  6 
/0 01 222 
/022122 

/ 0  0 0  24  0 

/0W5267 
Z000524 
/O01212 
/016700 
,-'0  0  0462 
/  0  2  0  2  6  7 
/ 0  0  0  4  6  4 
/ 10  0  4  0  6 
/  0 1  6  7  0 1 
,■'  0  0  0  4  5  2 
/ 016  7  0  2 

/ 0  0  0  4  5  0 

/162707 
/  0  0  0  7  0  0 
/  0 1 6  7  0 1 
/ 0  0  0  4  4  2 
/0167B2 
/  0  0  0  4  4  0 
/ 16  2  7  0  7 
/  0  0  0  7 1 4 

/  0  0  0  0  0  0 
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APPENDIX  E 
Program  Listing  of  Supplementary  Subprograms 

1.   Initialization  subprogram  (10174-10706) 

This  subprogram  enters  the  necessary  addresses  and 
constants  into  those  locations  noted  in  Appendix  A  by  * . 
The  computer  listing  below  is  the  initialization  subprogram. 


010174 
010176 

01020  0 

010202 
010204 
0 1 0  2  0  6 

010210 

010212 
010214 
010216 
010  2  2  0 
010  222 
010224 
0102  26 
010  230 
010232 
010234 
0 1 0  2  3  6 
010  240 
010242 
010  2  4  4 
010  24  6 
0 1 0  2  5  0 
0 10  2  5  2 
010254 
01025  6 
010260 
010  26  2 
010264 
010  2  6  6 


/010700 

,-•  0  1. 2  7  6  7 
/  0  0  0  0  0  2 
,'0034  74 
/  0 1 2  7  6  7 

/  0  0  0  0  0  3 

/ 0  0  3  4  7  0 
/ 012  7  6  7 

/  0  0  0  0  0  1 

/  0  £1  3  4  6  4 
,'012767 
/  0  0  0  0  2  0 
/  0  0  3  4  6  0 
/ 012  7  6  7 
/  0  0  0  0  0 1 
/ 0  0  3  4  5  6 
/ 012  7  6  7 
/  0  0  0  0  0  0 
/ 0  0  3  4  5  2 
/  0 1 0  7  0 1 
/0627EU 
/  0  0  3  5  1  4 
/ 01016  7 
/ 0  0  3  3  5  2 
/  0  6  2  7  0 1 

/  0  0  0  4  0  6 

/ 0101 6  7 
/ 0  0  2  3  4  4 
/  0  6  2  ?  0  0 
/ 0  0  3  4  0  2 


010270 

/  0  1  0  0  0  1 

010272 

J  0  6  2  7  0  0 

010274 

/  0  0  0  6  0  0 

010276 

/010021 

010300 

/&€■€■?  0  0 

010  30  2 

/ 0  0  3  4  0  4 

010304 

/  B 1 0  0  2 1 

610306 

/016737 

010310 

/  0  0  3  3  6  6 

010312 

/ 1  7  7  3  0  4 

010314 

/016737 

010316 

/  0  0  3  3  7  0 

010320 

/I  77306 

010322 

/013767 

010324 

/177304 

010326 

Z003420 

010330 

/  0  6  3  7  0  0 

010332 

/177304 

010334 

/010021 

010336 

/ 016  7  3  7 

010340 

Z003342 

010342 

/ 1  7  7  3  0  6 

010344 

/013767 

010346 

/ 17  7  3  0  4 

010350 

/  0  0  3  3  4  6 

0 10  3  5  2 

Z062700 

010354 

/177204 

0 1 0  3  5  6 

/  0  1  0  0  2 1 

010  36  0 

/  0  6  3  7  0  0 

0 1 0  3  6  2 

•177204 

90 


010364 
010366 
010370 
010372 
010374 
010376 

010  4  0  0 

010402 
0104EH 
010406 

010410 

010412 
010414 
010416 
010420 
010422 
010424 
010  4  2  6 
010430 
010432 
010434 
010  4  3  6 
010  4  4  0 
01044  2 
010444 
010446 
010450 
010  4  5  2 
010454 
010456 
010460 
010462 
010464 
010466 
010470 
010472 
010474 
O10  4  7  6 

01050O 

010  50  2 
010504 
010  5  0  6 

010  510 

010512 
010514 
010  516 
010  5  2  h 
0 1 0  5  2  2 
010  5  2  4 

010526 
O10  5  3  0 
010  5  3  2 
010  5  3  4 


•019021 

•  0  6  2  7  0  0 
•177304 

•  0  1  0  0  2 1 

/  0  0  5  0  0  2 

•005202 

•  0 1  6  7  0  3 
•002276 

•  0  Q  6  3  0  2 
•077302 
• 010  26  7 

•002210 

•  0  0  6  3  0  2 

•  0  6  0  2  0  0 

•  0  1  0  0  2  1 

•  0  0  6  2  0  2 

•  0  0  5  3  0  2 
•01026  7 

•  0  0  2  2  7  4 

•  0  0  5  2  8  2 

•  0 1  2  7  0  3 
/  0  0  0  0  0  2 

•  0  6  0  2  0  0 

•  0  1  0  0  2  1 

•077303 

•  0  0  6  2  0  2 
•010267 

•  0  0  2  2  5  4 
•005292 
•010267 
•002250 

•  0  0  5  2  0  2 

•  0  0  6  2  0  2 

•  010  2  6  7 

•  0  0  2  2  4  2 

•  012  7  2  7 
/  0  0  0  0  0  4 
•177304 
•016727 
/082176 
•177206 

•  0  1  2  7  6  7 

•  1  7  7  2  0  4 

•  0  0  2  2  0  6 

•  0 1 6  7  2  7 

•  0  0  2  1  6  4 
•177204 
•016727 

•  0  0  2 1 5  2 
•177206 
•012767 
/177204 

•  0  0  2  1 6  0 


010526 
010540 
010  5  4  2 
010544 
010546 
010550 
010552 
010554 
010  5  5  6 
010  5  6  0 
010562 
010564 
010566 

0105  70 
010572 
010574 
010576 

010600 

010602 
010604 
010  606 

010610 

010612 
010614 
01O616 
010620 
010  622 
010624 
0 10  6  2  6 

010  620 
010622 
010624 
010636 
010640 
010642 
010644 

0 106  4  6 
0 10  6  5  0 
010652 
010654 
0 1 0  6  5  6 
0 1 0  6  6  0 
0106  6  2 
010664 
010666 
010  6  7  0 
0 10  6  7  2 
010674 
010676 

010  7  0  0 

010  7  0  2 
010704 
010  7  0  6 


•  0 1 6  7  0  2 
•003142 

•  0  0  5  8  0  2 

•  0  0  5  2  0  2 

•  0  0  6  2  0  3 

•  0  7  7  2  0  2 

•  0  0  5  2  0  3 
•010267 

•  0  0  2 1  5  6 

•  0  0  5  2  0  2 

•  0  0  6  2  0  2 
•010267 

•  0  0  210  6 
•016702 

•  0  0  2 1 0  6 

•  0  0  5  2  0  2 

•  0  6  6  7  0  2 
•002104 
•010267 
•002122 

•  0 1  6  7  0  2 

•  0  0  2  0  7  4 

•  0  0  5  2  0  2 

•  0  0  5  0  0  3 

•  0  0  5  2  0  2 

•  0  0  6  2  0  2 

•  0  7  7  2  0  2 
•01026  7 

•  0  0  211 2 

•016702 

•  0  0  2  7  6  4 
•166702 
•002756 
•010267 

•  0  0  2 1 0  0 

•  0 1 6  7  0  2 
Z062752 

•  16  6  7  0  2 

•  0  0  2  7  4  2 

•  0 1 0  2  6  7 

•  0  0  2  0  6  6 

•  0  1  0  7  0  0 

•  0  6  2  7  0  0 

•  0  0 12  6  4 

•  0 1 0  0  6  7 
/  0  0  2  2  0  4 
•006267 

•  0  0  2  0 1 2 
• 0  0  6  2  6  7 

•  0  0  2  0  1 6 

•  0  0  6  2  6  7 
/ 0  0  2  0  4  2 

/  0  0  0  0  0  0 
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2.   Message  Input  Subprogram  (10000-10144) 

This  subprogram  was  used  to  indicate  that  an  alphanu- 
meric symbol  typed  at  the  Keyboard  was  entered  into  core 
(echo) .  The  symbols  are  stored  in  core  in  an  ASC  II 
(7-bit)  code  representation.   The  Keyboard  symbol,  @  , 
is  used  to  terminate  message  entry  and  the  number  of 
8-bit  computer  bytes  used  is  stored  in  location  (13710 
Appendix  A) .   The  following  computer  printout  is  a  listing 
of  the  ASC  II  message  input  subprogram. 


016000  7010700  610864  7600012 

010002  /0G2700  010066  7105727 

010004  7004376  010076  7177564 

010006  7005062  010072  7160275 

010010  7105737  010074  7112737 

010012/177560  010076  .'000200 

010014  /l 002  7  5  010100  7177566 

010016  7113  710  01010  2  /0 7 7 10 7 

010020  7177562  016104  7105737 

010022  /l 22710  010106  7177564 

010024  7000300  010110  7100  375 

010026  7001435  010112  7112  737 

010030  7105737  010114  7000212 

010032  7177564  010116  7177566 

010034  71E-0375  010120  700  0  7  33 

010036  7112037  010122  7012703 

01004  0  717  7  566  010124  7606010 

01004  2  7  005202  010126  7105626 
61004  4  7123  727  010136  7065262 
010  64  6  7177562  01013  2  70  7  7  30  2 
016  650  70  60  215  O10134  76  66262 

01005  2  760135  6  016136  70  6  6  302 
0100  5  4  7112  740  016146  7610267 

01005  6  7006  24  6  01014  2  7  0  03544 

01006  0  710  5  726  610144  706  6666 
010062  7612761  * 
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3.   Message  Output  Subprogram  (13200-13364) 

This  subprogram  is  basically  the  same  as  that  in 
section  2  of  this  appendix.   However,  besides  typing  what 
is  in  core  starting  at  a  location  specified  at  (13602) , 
the  number  of  errors  (bit  differences)  are  determined 
between  the  input  message  and  the  decoded  out  message. 
This  message  error  value  is  stored  at  (13672)  (Appendix  A) 
The  following  computer  printout  is  a  listing  of  the  ASC  II 
(7-bit)  code  output  subprogram. 


012200 
012202 
012204 
01220  6 
012210 
012212 
012214 
012216 
01222  0 
012  22  2 
012224 
012226 
012220 
012222 
0 1 2  2  2  4 
012226 
012240 
012242 
012244 
0122  4  6 
CU2250 
012252 
012254 
012256 
0122  6  0 
012  26  2 
0 1 2  2  6  4 
012266 
012270 
012272 


/016700 
/ 0  0  0  2  7  6 
/  012  7  01 

Z000040 

•111027 
,'177566 
•122720 

•  0  0  0  0  0  0 
• 0  014  2  5 

• 10  5  7  2  7 

•  17  7  5  6  4 
• 10  0  2  7  5 

•  0  0  5  2  01 
• 1 0  0  2  6  6 
•122710 

•  0  0  0  2  4  0 
/  8  0 1 2  6  2 

.'A     A     *"i  "?  -.-  -> 

f  i  1  <i  r  S  I 

• 0  0  0  015 

•17  7  5  6  6 
• 012  7  0  2 

•  0  0  0  012 
/ 10  5  7  2  7 
• 1  7  7  5  6  4 
/ 10  0  2  7  5 
•185037 
/ 177566 

•  0  7  7  2  0  6 
•105727 

•  1 7  7  5  6  4 


012274 

• 1 0  0  2  7  5 

012276 

• 112  7  2  7 

61220  0 

/000012 

013202 

/ 1  7  7  5  6  6 

012204 

/ 1 0  5  7  2  7 

013306 

/ 1 7  7  5  6  4 

013310 

/100275 

013212 

/  0  0  0  7  2  4 

012214 

Z000  24  0 

013316 

/ 0  00 240 

013320 

/  0 1 6  7  0  0 

013  2  2  2 

/ 0  0  0  2  5  4 

012224 

/ 0 1 6  7  01 

012326 

/ 0  0  0  2  5  2 

013  32  0 

/  0  0  5  0  0  2 

81iii2 

/011002 

012224 

/ 0  7  4  211 

012226 

/  0 1 2  7  0  4 

01334  0 

/  0  0  0  0  2  0 

013342 

/ 0  0  6  211 

013244 

/  0  0  5  5  0  2 

012  2  4  6 

/  0  7  7  4  0  2 

0 1 2  2  5  0 

/  0  2  2  0  2 1 

012252 

/  0  2  0  0  6  7 

012  2  5  4 

/  0  0  0  2  2  4 

012256 

/001265 

0 1 2  2  6  0 

/ 010  2  6  7 

012262 

/  0  0  0  2  0  6 

012  2  6  4 

•000000 

* 
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4.   Analysis  Subprograms  (13250-13544) 

From  address  (13250)  to  address  (13376)  a  short  sub- 
program is  listed  that  is  used  to  step  through  well  defined 
encoders  of  rate  1/2,  1/3,  and  1/4.   In  the  last  addresses 
(13400-13544)  a  program  is  listed  which  determines  the 
number  of  error  bits  in  16-quantized  bits  and  stores  the 
results  so  that  a  distribution  of  errors  may  be  plotted 
as  in  figures  17-23.   (See  next  page.) 
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013250 
01225  2 
012254 
012256 
012260 
012262 
012264 
012266 
012270 
B13272 
012274 
012276 
0122  0  0 
012202 
012204 
012  20  6 
012210 
012212 
012214 
012216 
012220 
012222 
012224 
012226 
012220 
012222 
012224 
012  226 
012240 
0122  4  2 
012  24  4 
012246 
012250 
012252 
012254 
012  25  6 
012  26  0 
B12262 
012264 

012  2  6  6 
012270 

013  2  7  2 
012274 
012  2  7  6 


•005267 

/  0  0  0  5 1 2 
/ 0  2  6  7  2  7 

• 0  0  0  5  0  6 
/000010 
/0010  20 

•012767 

•  0  0  0  0  0  4 

• 0  0  0  4  7  4 
/  0  0  5  2  6  7 

•  0  O  0  4  6  6 
/ 026 727 

•  0  0  3  4  6  2 

•000010 
/ 0  01017 

•  0 1 2  7  6  7 

/  0  0  0  0  0  4 
• 0  0  0  4  5  0 

•005267 

•  0  0  0  4  4  2 

•  0  2  6  7  2  7 
/  0  0  O  4  2  6 

•  0  0  0  0  2  0 

• 0  0 10  0  6 

/012767 

•  0  0  0  0  1  0 

• 0  00 4  24 
/  0  6  2  7  6  7 

•  0  0  0  0  0  2 

•  0  0  0  414 
/016777 
•000222 

•  0  0O416 

•  0  6  2  7  6  7 

/  0  0  0  0  0  2 
/ 0  0  6  410 

•  0  2  6  7  2  7 

•  0  0  0  2  7  4 
,-'  0  0  ij  0  2 1 
/  0  014  0  2 

•  0  0  012  7 

/01000  0 

/  0  0  9 1 2  7 
/  0  0 1 1  7  2 


0134  00 
012402 
012404 
012406 
012410 
012412 
012414 
012416 
612420 
012422 
012424 
012426 
0124  20 
012422 
012424 
0124  2  6 
012440 
012442 
012444 
012446 
012450 
012452 
012454 
012456 
012460 
012462 
012464 
012466 
GU2470 
012472 
012474 
013476 
012500 
012502 
012504 
012506 
O12510 
012  512 
012514 
012516 
012520 
012522 
012  524 
01252  6 
012520 
012522 
012524 
012  5  26 
012  54  0 
01254  2 
012  5  4  4 


•016700 
•000202 
•016781 
•000226 
/  0  0  6  2  0 1 

•  0  O  0  2  4  0 
/  012  7  6  7 

/  0  0  0  0  2  0 
•000214 

•012784 

•014108 

/012705 

•  0  0  0  0  2  0 
,'0  0  0  24  0 

•080240 
•800240 

•000240 

•  0  0  5  0  0  6 

•016703 

•000224 

•  0  0  5  0  0  2 
•006282 

•  0  0  6  2 1 0 
•005582 
•005367 
•000152 

•  0  0 1 4  0  7 
•077387 
•020201 

•  0  0  2  4  81 

•  0  0  5  2  0  6 
•005205 

•  0  0 1 4 1 4 

•  O  &&?6 0 
•0  05  720 
•0  26  70  0 

•  0  0  0  0  7  6 

•  0  0  2  0  0  2 

,'  0  0  0  0  0  0 

•  0  0  0  2  4  0 

•  0  0  0  2  4  0 

•012767 

•  0  0  O  0  2  0 

•  0  0  010  6 

•  0  0  0  7  5  6 

•  0  0  6  2  0  6 

•  0  6  0  6  0  4 

•  0  0  5  214 

•  0  0  0  2  4  0 

•000727 

•  0  0  O  O  0  0 
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APPENDIX  F 
SAMPLE  RUN 

THE  COMPUTER  LISTING  BELOW  THIS  PARAGRAPH  IS  A 
DEMONSTRATION  OF  THE  PROCEDURE  REQUIRED  TO  USE  THE 
ENTIRE  PROGRAM  PROPERLV.  THE  PROGRAM  IS  STORED  ON 
A  DISK  UNDER  THE  NAME  VSAC.  SRV.  THE  CONVOLUTI ONAL 
CODE  USED  IS  A  RATE  1/2- (K=3>  CODE.  OTHER  INPUTS 
ARE  Q  =  l,  DCL=16  <C20  IN  BASE  8),  FIND  THE  GENERATOR 
SEQUENCES  OF  3  AND  4  (BASE  8). 


.  GET  VSAC.  SRV 


.  STAR 

ODT 

* 

102O 
1020 

1021 
1022 
1023 
1023 


137b 
1376 


T  1172 

V01-01 

0/000062 
fc"  /  0  0  0  0  0  3 
4  /  0  0  0  U  0  1 
2/0000  20 
0  /  0  0  0  0  0  1 

:6/000O01 
:6/100004 

0  /  0  0  0  0  0  3 
2/000004 


"C 

START  10000 
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The  paragraph  below  is  the  input  message  for  the  sample 


run. 


THIS  IS  R  SRMPLE  RUN  OF  ft  THESIS  COMPUTER 
PROGRAM  USED  TO  EVALUATE  THE  PERFORMANCE  OF  RATE 
i/N  CONVOLUTIONAL  CODES,  OVER  A  SIMULATED  NOISV 
CHANNEL,  WITH  A  VITERBI  DECODER.    THE  NOISE  PARAMETERS 
FOR  THIS  RUN  <CORE  LOCATIONS  10230  AND  10236)  PRE 
THE  SftME  AS  WERE  USED  TO  OBTAIN  THE  DISTRIBUTION 
IN  FIGURE  20  OF  THE  THESIS  TEXT.    THIS  MESSAGE  WILL 
NOW  BE  TERMINATED  BV  THE  'AT'  CHARACTER  ON  THE  PANEL. 


When  the  message  above  is  encoded  (as  given  above)  and  the 

noise  is  added,  the  resulting  decoded  message  is: 

THIS  IS! ft"SGMPLE  RUN  OF  E  THESIS 
COMPUTR  PROGBAM  USED  \0  EVALUATE 
THE  ERFNRURNCEBDF  RRTE  1/N  CONVOLUTIONAL ! BCDES, 
OVMR  RPIMULBTEE  NOISV  PRNNEL   ITHCG 
VITERBI  DECOLDP.    UHE  NOISE  PARAMETERS 
FOR  TPIS  UN  (COU  LOCATIONS  18228 
AND  102261  RRU  $E    SRUR  RS  TEVE  USED" TO 
OBTAIN  THE  DISTRIBUTION  INBFJURE 
21  OF  THE  THE  IP  TEXD.    THIS  MESSRGD 
WILL  NIW  BE  TERMINATED  BV'THM  'AX' 
CHBRCTER  ON  THM  BBNEB. 

When  no  coding  is  applied  to  the  input  message  and  the 

same  noise  is  added  the  received  message  is: 

IB IL  fl I? 0 [<  I  N B  J R G 8 Elfl  J B M E  B  V U N  "  M D 8  0.  (.  T  I  C K C  8 C 0 M X U 
D  R  R  0  A  M  <  U  R  M  B  *  D  _  <  M  T  C  M  1  fl  U  E  D  B  R  8  w  D  R  F  0  E  R  F  K  D  "  F 
S  K  U I  "  A  !  R  H  V  R  V  R  H  "  E  G  C  0  E  R  >  8  B  \  B  D  !  N  N  R  S  D  "  T  C  R  Q I 
UUGRSONKV  <r- •  * 

PJIS  ZUO  !  HCORESLBITIOJC  10  :  7"ED"943i  V9* 
R  &  E  B  I  0  U  $  I  B  f  G  R  M  8  E  C  U  D  B  T  0  8  M  B  V  E I  N  !  D  U  D  S  T  P  H  C  W  U  I 
VDffEBJEBQIRU!SQT:L:fNOG<E  TGR 
BJRIKPRR"_NDJM! KNEL.  @ 

*  Improper  carriage  return  received. 
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